1/2 


OD-HlSt  279  DEVELOPMENT  OF  R  MODEL  FORRTHE  HERT  RELERSE  RRTE  OF 
UUOD  R  STATUS  REPORT(U)  NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON  DC  CENTER  FOR  FIRE  RE. .  M  J  PARKER  HRV  89 
UNCLASSIFIED  NBSIR-89-2163  EHH-E-1239  F/G  11/12 


jLvv-v  v  %  v.vw;. 


NBSIR  85*3163 


Development  of  A  Model  for  the 
Heat  Release  Rate  of  Wood  - 
A  Status  Report 


William  J.  Parker 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standarda 
National  Engineering  Laboratory 
Cantor  for  Fire  Research 
Gaithersburg,  MD  20899 


May  1986 


DT1C 

!  ELECT E 
.  OCT  1 6  B85 


OTIC  file  coey 


Sponsored  in  part  by: 

Federal  Emergency  Management  Agency 


This  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited 


85  10  15  04  5 


Nesm  88-3163 


DEVELOPMENT  OF  A  MODEL  FOR  THE 
HEAT  RELEASE  RATE  OF  WOOD  - 
A  STATUS  REPORT 


William  J.  Parfcar 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standards 
National  Engineering  Laboratory 
Center  for  Fire  Research 
Gaithersburg,  MD  20699 


May  1985 


Sponsored  in  part  by: 

Federal  Emergency  Management  Agency 


Acce^ 

ion  For 

NT!  3 

CRA&I 

t,  : 

TAB 

U  u  .■ 

iOLJ,  .COC! 

J  i.i.Cution 


By 

n  .  lion/ 


f  *  nubility  Codes 


I  Av«.'il  andjor. 
i  Sise  oaf 

:  I 


U.S.  DEPARTMENT  OF  COMMERCE,  Malcolm  Baldrige,  Secretary 

NATIONAL  BUREAU  OF  STANDARDS,  Ernest  Ambler.  Director 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (Whan  Data  Cntarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  1 

BEFORE  COMPLETING  FORM  I 

nEsaaa  **  n-i.'rni 

«.  TITLE  (mt4  Submit) 

DEVELOPMENT  OF  A  MODEL  FOR  THE  EE AT 

RELEASE  RATE  OF  WOOD  -  A  STATUS  REPORT 

3.  type  of  report  a  period  covereo 

Final 

«.  PERFORMING  ORG*  REPORT  NUMBER 

NBSIR  85-3183 

7.  AUTHORS 

William  J.  Parker 

i.  CONTRACT  OR  GRANT  NUMBER!-*; 

EMW-E-1239 

S.  PERFORMING  organization  name  ano  AOORESS 

Center  for  Fire  Research 

U.S.  National  Bureau  of  Standards 
Gaithersburg,  MD  20899 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

2531  G 

It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Federal  Emergency  Management  Agency 
Washington,  D.C.  20U72 

12.  report  date 

May  1985 

11.  NUMBER  OF  PAGES 

106 

<4.  MONITORING  AGENCY  name  A  AOORESSfi!  illlatanl  tram  Cantralllna  0111  ca) 

•  S.  SECURITY  CLASS,  (ol  title  report) 

Unclassified 

ISa.  OECLASSI F| CATION/  OOWNGRAOING 

SCH SOULE 

IS.  DISTRIBUTION  STATEMENT  (al  tbit  Kepart) 

Approved  for  public  release:  Distribution  Unlimited 

17.  DISTRIBUTION  STATEMENT  (al  the  abaUaet  antarai  In  Black  20.  II  illlatanl  tram  Bap an) 

Approved  for  public  release:  Distribution  Unlimited 

is.  supplementary  notes 

1*.  KEY  WORDS  (Continue  on  rooeroe  «ido  If  nocoooory  Identify  by  block  number) 

heat  release  rate,  wood,  pyrolysis,  heat  of  combustion, 
thermal  conductivity,  computer  models 

20}  ABSTRACT  rramttaua  am  not  lip  U  nmamaatf  ami.  Idanllfy  by  black  number) 

This  report  describes  the  status  of  the  development  of  a  method  foi 
predicting  the  heat  release  rate  of  wood  for  different  thicknesses, 
moisture  contents,  and  exposure  conditions.  A  computer  program  has 
been  set  up  on  a  microcomputer.  Experimental  techniques  have  been 
devised  to  obtain  the  input  data  required  by  the  model.  These  in¬ 
clude  (1)  the  thermal  conductivity  as  a  function  of  temperature  anc 
percent  loss,  (2)  the  kinetic  constants  needed  to  describe  the  mass 

00  waTts  M73  COITION  OF  I  NOV  SS  IS  oosocete 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  PAGE  D*f*  Cntattt)) 


UNCLASSIFIED _ 

IKCmWTY  OMBHCATjOII  OF  THIS  FAOKlWNaH  Qmm  Mmtan4) 


(Block  20) 

,loss  rate,  (3)  the  heat  of  combustion  of  the  volatile  pyrolysis 
products,  and  (k)  the  contraction  factors  due  to  charring.  Suffi¬ 
cient  data  on  these  parameters  were  taken  to  exercise  the  model. 
Heat  release  rates  and  effective  heats  of  combustion  were  measured 
as  a  function  of  external  radiant  flux  on  12.5  nun.  thick  dry  verti¬ 
cal  specimens  of  Douglas  fir  particle  board.  The  calculated  and 
measured  peak  heat  release  rate  curves  are  similar  in  shape  and 
amplitude  but  differ  significantly  in  time  scale.  This  may  Be  due 
to  the  lack  of  thermal  conductivity  data  on  the  char  in  the  high 
temperature  range.  There  is  very  good  agreement  Between  the  cal¬ 
culated  and  measured  effective  heats  of  combustion.  The  initial 
results  with  the  model  are  promising. 


UNCLASSIFIED 


MCUftITY  CLASSIFICATION  OF  THIS  FAOKPFhOT  Data  fntara*) 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES .  v 

LIST  OF  TABLES  .  viii 

NOMENCLATURE .  lx 

Abstract  . . . . . . .  1 

1.  INTRODUCTION . . .  1 

2.  HEAT  RELEASE  RATE  MODEL .  5 

3.  THERMOCHEMICAL  PROPERTIES  .  15 

3. 1  Pyrolyzer  .  15 

3.2  Heat  of  Combustion  of  Pyrolysis  Vapor . 17 

3.3  Calculation  of  Kinetic  Constants  from  Mass  Loss  Data  at 

Constant  Temperature . . . 18 

3.4  Contraction  Factors . . .  20 

3.5  Heat  of  Pyrolysis . 20 

3.6  Parameters  Dealing  with  Moisture  .  21 

4.  THERMOPHYSICAL  PROPERTIES  .  21 

4. 1  Heat  Capacity  .  21 

4.2  Thermal  Conductivity  Data  from  Literature  .  22 

4.3  Thermal  Conductivity  Measurements  . .  24 

5.  CALCULATIONS  OF  HEAT  RELEASE  RATE  USING  THE  MODEL . 28 

6.  HEAT  RELEASE  RATE  MEASUREMENTS  IN  THE  CONE  CALORIMETER .  30 

7.  DISCUSSION .  32 

8.  SUMMARY  .  40 

9.  FUTURE  PLANS .  41 

10.  ACKNOWLEDGEMENTS .  42 

11.  REFERENCES .  43 

APPENDIX  A.  CALCULATION  OF  THE  DENSITY  OF  THE  GAS  IN  THE  EXHAUST  DUCT  ..  52 

APPENDIX  B.  CALCULATION  OF  THE  FLOW  OF  NITROGEN  THROUGH  THE  EXHAUST 

DUCT  . .  55 

APPENDIX  C.  CALCULATION  OF  THE  CARBON,  HYDROGEN,  AND  OXYGEN  RELEASE 

RATES  FROM  BURNING  WOOD . 57 


ill 


LIST  OF  FIGURES 

Pag 

Figure  1.  Boundary  Conditions  and  Subdivision  of  Burning  Vertical  Slab  6 


Figure  2.  Flow  Chart  for  Computer  Program .  6 

Figure  3.  Pyrolyzer . 6 

Figure  4.  Closure  Slide  . . 6 

Figure  5.  Heat  of  Combustion  of  Pyrolysis  Vapor  for  Douglas  Fir 

Particle  Board . . . . . . .  6 

Figure  6.  Arrhenius  Plot  for  Douglas  Fir  Particle  Board  .  7' 

Figure  7.  Contraction  Factors  for  Southern  Pine  . .  7 

Figure  8.  Contraction  Factors  for  Douglas  Fir  Particle  Board . .  7 

Figure  9.  Experimental  Setup  for  Determining  Thermal  Conductivity  .  7 

Figure  10.  Thermal  Conductivity  as  a  Function  of  Temperature  for 

Charred  and  Uhcharred  Specimens  of  Southern  Pine  .  7 

Figure  11.  Thermal  Conductivity  as  a  Function  of  Temperature  for 


Charred  and  Uncharred  Specimens  of  Douglas  Fir  Particle  Board  7 

Figure  12.  Calculated  Front  and  Rear  Surface  Temperatures  for  Dry  12.7  mm 
Specimen  of  Douglas  Fir  Particle  Board  Exposed  at  an  External 


Radiant  Flux  of  50  kW/m2  . .  7 

Figure  13.  Calculated  Char  Depth  for  Dry  12.7  mm  Specimen  of  Douglas  Fir 
Particle  Board  Exposed  at  an  External  Radiant  Flux 
of  50  kW/m2  .  7 

Figure  14.  Calculated  Mass  Loss  Rate  for  Dry  12. 7  mm  Specimen  of  Douglas 
Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 
of  50  kW/m2 .  7 

Figure  15  Calculated  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of 
Douglas  Fir  Particle  Board  Exposed  at  an  External  Radiant 
Flux  of  50  kW/m2  . . .  7 

Figure  16.  Calculated  Effective  Heat  of  Combustion  for  Dry  12.7  mm 

Specimen  of  Douglas  Fir  Particle  Board  Exposed  at  an  External 
Radiant  Flux  of  50  kW/m2  .  8< 

Figure  17.  Calculated  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 
Fir  Particle  Board  with  an  Assumed  Heat  of  Pyrolysis  of  400  kJ/kg 
Exposed  at  an  External  Radiant  Flux  of  50  kW/m2  .  8 


List  of  Figures  (continued) 


Figure  18. 


Figure  19. 


Figure  20. 


Figure  21. 


Figure  22. 


Figure  23. 


Figure  24. 


Figure  25. 


Figure  26. 


Figure  27. 


Figure  28. 


Figure  29. 


1*S£ 

Calculated  Heat  Release  Rate  for  IX  Moisture  Content  12.7  on 
Specimen  of  Douglas  Fir  Particle  Board  Exposed  at  an  External 
Radiant  Flux  of  50  kW/m2 .  82 

Calculated  Heat  Release  Rate  for  Dry  6  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  50  kW/m2  .  83 

Calculated  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  25  kW/m2  . 84 

Calculated  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 
Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux  of 
75  kW/m2  .  85 

Calculated  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux  of 

100  kW/m2  .  86 

Measured  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  25  kW/m2  .  87 

Measured  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  50  kW/m2  .  88 

Measured  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  75  kW/m2  .  89 

Measured  Heat  Release  Rate  for  Dry  12.7  mm  Specimen  of  Douglas 

Fir  Particle  Board  Exposed  at  an  External  Radiant  Flux 

of  100  kW/m2  .  90 

Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Fir  Particle  Board  Exposed  at  an  External  Radiant 
Flux  of  25  kW/m2  .  91 

Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Fir  Particle  Board  Exposed  at  an  External  Radiant 
Flux  of  50  kW/m2  .  92 

Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Fir  Particle  Board  Exposed  at  an  External  Radiant 
Flux  of  75  kW/m2  .  93 


List  of  Figures  (continued) 


Figure  30.  Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Fir  Particle  Board  Exposed  at  External  Radiant 
Flux  of  100  kW/2  . 


Figure  31. 


Measured  Heat  Release  Rate  Determined  by  the  Three  Analyzer 
Method  for  Dry  12.7  mm  Specimen  of  Douglas  Fir  Particle 
Board  Exposed  at  an  External  Radiant  Flux  of  50  lcW/m2  ...... 


Figure  32.  Measured  Effective  Heat  of  Combustion  Determined  by  the  Three 
Analyzer  Method  for  Dry  12.7  mm  Specimen  of  Douglas  Fir 
Particle  Board  Exposed  at  an  External  Radiant  Flux 
of  50  ktf/m2  . . 


LIST  OF  TABLES 


Thermal  Conductivity  Data  on  Flberboard  and  Particle  Board 
from  the  Literature. . . . 

Thermal  Conductivity  of  Oven  Dry  Southern  Pine  as  a  Function 
of  Temperature . . . . . 

Thermal  Conductivity  of  Charred  Southern  Pine  as  a  Function 
of  Temperature . . . 

Thermal  Conductivity  of  Oven  Dry  Douglas  Fir  Particle  Board  as 
a  Function  of  Temperature . . . . . . 

Thermal  Conductivity  of  Charred  Douglas  Fir  Particle  Board  as 
a  Function  of  Temperature . . . 

Input  Data  for  Running  the  Model  for  Douglas  Fir  Particle 


NOMENCLATURE 


frequency  factor  (s“*) 
mass  transfer  number 
heat  capacity  (IcJ/kgK) 
activation  energy  (kJ/mol) 

heat  of  combustion  of  volatiles  (kJ/kg)  (positive  when  heat  is 
evolved) 

convective  heat  transfer  coefficient  (ktf/m^K) 

heat  of  pyrolysis  (kJ/kg)  (positive  when  heat  is  absorbed) 

thermal  conductivity  (kW/m  •  K) 

contraction  factor 

effective  heat  of  gasification  (kJ/kg)  (positive) 
latent  heat  of  vaporization  of  water  (kJ/kg)  (positive) 
mass  (kg) 

Z  moisture  content 
number  of  components 

number  of  subdivisions  of  the  specimen 

molar  flow  (mol/s) 

grouping  defined  by  equation  (2-14) 

heat  (kJ)  (positive) 

external  radiant  flux  (kW/m^) 

rate  of  heat  release  (kW) 

gas  constant  (8.314  kJ/mol  •  K) 

time  (s) 

temperature  (°C) 
volume  flow  (m^/s) 

distance  up  from  bottom  of  specimen  (m) 


volume  fraction 


mass  retention  fraction 

radiative  absorptivity  and  emissivlty 

heat  of  combustion  (kJ/kg)  (positive  when  heat  is  evolved) 
pressure  drop  across  orifice  plate  (Pa) 
time  step  (s) 

thickness  of  each  slice  (m) 
defined  by  equation  (2-20) 
density  (kg/nr*) 

Stef an-Boltzmann  constant  (5.67  x  10“**  kW/m^  •  K^) 
ts 

convection 

carbon 

thickness  of  specimen  (m) 
external  or  exhaust 
flame 
gas 

hydrogen 

designates  a  particular  slice 

designates  a  particular  wood  component 

designates  the  slice  bounded  by  the  rear  surface 

ambient  or  original 

oxygen 

radiation 

release 

front  surface  of  specimen 


volatiles 


w 


water 


x  normal  to  surface 

y  parallel  to  surface  normal  to  grain 

z  parallel  to  surface  along  grain 

Superscripts 
A  analyzer 

A'  analyzer  with  CO2  trapped  out  of  it 

S  exhaust  duct 

•  per  unit  time 


per  unit  area 


DEVELOPMENT  OF  A  MODEL  FOR  THE  HEAT  RELEASE  RATE  OF  WOOD  -  A  STATUS  REPORT 


William  J.  Parker 

Abstract 

This  report  describes  the  status  of  the  development  of  a  method  for 
predicting  the  heat  release  rate  of  wood  for  different  thicknesses,  moisture 
contents,  and  exposure  conditions.  A  computer  model  has  been  set  up  on  a 
microcomputer.  Experimental  techniques  have  been  devised  to  obtain  the  input 
data  required  by  the  model.  These  include  (1)  the  thermal  conductivity  as  a 
function  of  temperature  and  percent  mass  loss,  (2)  the  kinetic  constants 
needed  to  describe  the  mass  loss  rate,  (3)  the  heat  of  combustion  of  the 
volatile  pyrolysis  products,  and  (4)  the  contraction  factors  due  to  charring. 
Sufficient  data  on  these  parameters  were  taken  to  exercise  the  model.  Heat 
release  rates  and  effective  heats  of  combustion  were  calculated  and  measured 
as  a  function  of  external  radiant  flux  on  12.5  mm  thick  dry  vertical  specimens 
of  Douglas  fir  particle  board.  The  Initial  results  with  the  model  are 
promising. 

Key  words:  heat  release  rate,  wood,  pyrolysis,  heat  of  combustion, 
thermal  conductivity. 


1.  INTRODUCTION 


The  heat  release  rate  is  one  of  the  most  important  fire  properties  of  a 
material.  Whether  a  room  fire  will  attain  flashover  depends  on  the  total  heat 
release  rate  of  all  the  materials  in  the  room  as  well  as  its  .ventilation  and 
surface  heat  losses  [1,2].  The  rate  of  flame  spread  in  the  direction  of  air 


flew  Is  strongly  dependent  on  the  rate  of  heat  release.  In  order  to  model  the 
fire  growth  In  a  room  it  is  essential  to  have  data  on  the  heat  release  rate  of 
all  of  the  materials  involved  as  a  function  of  the  net  heat  transfer  through 
their  front  surface. 

A  number  of  calorimeters  have  been  developed  to  measure  the  heat  release 
rate  per  unit  area  as  a  function  of  the  external  radiant  flux  [3.7].  The 
corresponding  net  heat  transfer  through  the  flaming  surface  as  required  by 
room  fire  growth  models  is  not  measured  in  any  of  these  calorimeters.  The 
heat  release  rate  curves'  as  a  function  of  time  at  a  particular  external 
radiant  flux  are  in  general  different  for  different  calorimeters  because  of 
the  different  conditions  under  which  the  data  are  obtained.  The  variables 
include  orientation,  specimen  size,  heat  transfer  conditions  at  the  rear  face, 
piloting  arrangements,  etc. 

An  approach  which  has  been  used  with  apparent  success  for  determining  the 
heat  release  rate  of  polymethylmethacrylate  (EMMA)  is  to  define  an  effective 
heat  of  gasification  as 

l  -  r„.e/*-  a-i) 

where  Q"net  Is  the  net  heat  transfer  through  the  surface  and  Si"  is  the  outward 
mass  flux.  The  mass  loss  rate  used  in  equation  (1-1)  is  measured  under  non¬ 
flaming  conditions  thus  avoiding  the  uncertainties  in  measuring  the  flame  heat 
transfer.  The  net  heat  transfer  to  the  specimen  is  equal  to  the  absorbed 
external  radiant  flux  minus  the  reradlatlon  and  convective  heat  losses  from 
the  front  face.  Tewarson  [8]  has  cataloged  values  of  L  and  the  effective  heat 
of  combustion,  (AH)  ,  for  many  materials. 


A  room  fire  growth  aodel  can  calculate  the  net  heat  transfer  rate  through 
the  surface  and  use  equation  (1-1)  along  with  the  experimentally  measured 
value  of  L  to  determine  the  mass  loss  rate.  The  calculated  heat  release  rate 
is  singly  equal  to 

«r.l  •  **  <*>*£•  <‘-2) 

Unfortunately  both  and  L  vary  throughout  the  burning  period  for  char 

forming  materials. 

When  L  reaches  a  steady  value  it  is  an  order  of  magnitude  greater  than 
the  heat  of  pyrolysis  of  cellulose  [9],  which  is  the  major  component  of  wood. 
This  suggests  that  L  in  this  case  is  dominated  by  heat  transfer  and  thus 
depends  on  density,  thickness,  etc.  Therefore,  the  effective  heat  of  gasifi¬ 
cation  is  not  a  useful  concept  for  wood. 

A  better  approach  for  char  forming  materials  is  to  develop  a  model  which 
predicts  the  mass  loss  rate  as  a  function  of  time  based  on  their  thermo¬ 
chemical  and  thermophysical  properties.  Two  general  types  of  mass  loss  rate 
models  have  been  developed.  Kanury  [10]  and  others  have  used  a  pyrolysis 
temperature  model  which  can  be  solved  analytically.  It  is  assumed  that  a  thin 
pyrolysis  front  moves  through  the  material.  Virgin  wood  exists  at  all  temper¬ 
atures  below  the  pyrolysis  temperature  and  the  final  char  exists  above  it. 

This  is  satisfactory  for  high  heat  fluxes  and  long  times. 

However,  for  fire  buildup  in  rooms  we  have  to  deal  with  the  whole  range 
of  fluxes  and  are  particularly  interested  in  early  times.  This  requires  the 
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use  of  kinetic  models  such  as  that  originated  by  Rung  [11]  and  inproved  upon 
by  Tamanini  [12]  and  Atreya  [13].  The  kinetic  models  allow  for  pyrolysis  to 
occur  at  any  temperature  at  a  rate  determined  by  an  appropriate  Arrhenius 
expression.  Again  there  are  a  number  of  other  models  but  it  Is  the  one 
refined  by  Atreya  which  will  be  built  upon  in  this  report. 

It  is  also  important  to  allow  for  moisture  evaporation.  The  temperature 
at  any  point  in  the  solid  is  determined  by  solving  the  energy  equation  numeri¬ 
cally. 

Even  the  kinetic  mass  loss  rate  models  are  not  without  problems.  Wood 
decomposes  along  two  principal  pathways.  In  one  case  all  of  the  gasified 
material  escapes;  it  will  generally  decompose  into  a  variety  of  products 
before  it  passes  through  the  surface.  It  will  have  an  overall  heat  of  combus¬ 
tion  similar  to  that  determined  with  an  oxygen  bomb  calorimeter.  Along  the 
other  pathway,  char  is  formed  with  the  release  of  water  and  other  volatiles 

having  low  heats  of  combustion.  The  effective  heat  of  combustion,  (aH)  ,,, 

erf 

will  be  determined  by  the  ratio  of  the  quantities  following  these  two  routes. 
This  ratio  changes  during  the  course  of  the  test.  The  mass  loss  rate  model 
cannot  accurately  predict  the  heat  release  rate  unless  the  change  in  effective 
heat  of  combustion  during  the  burning  period  is  known. 

In  this  report  a  model  is  presented  which  calculates  the  heat  release 
rate  directly.  Methods  of  providing  the  thermo chemical  and  thermophysical 
properties  required  as  input  data  for  the  model  are  described.  Some  of  these 
data  are  presented  for  Douglas  fir  particle  board  and  the  calculated  and 
measured  heat  release  rates  are  compared.  Douglas  fir  particle  board  was 


chosen  because  of  its  use  in  a  number  of  other  fire  research  projects.  The 
material  used  was  supplied  by  the  Weyerhaeuser  Company.  It  has  an  oven  dry 


density  of  709  kg/m^.  It  is  12.7  mm  thick  and  uses  urea  formaldehyde  (62  of 
the  mass  of  the  original  wood)  as  the  adhesive.  In  addition  there  is  1/4Z 
add-on  each  of  urea  and  a  paraffin-type  wax.  While  the  binder  may  have  some 
effect  on  the  chemical  and  physical  properties  of  the  particle  board,  these 
resultant  properties  are  measured  on  the  project. 

2.  HEAT  RELEASE  RATE  MODEL 


The  heat  release  rate  model  described  in  this  report  (1)  breaks  the 
specimen  down  into  thin  slices  of  equal  thickness  parallel  to  the  surface  as 
seen  in  figure  1,  (2)  calculates  the  teiq>erature  of  each  slice  using  the 
energy  equation,  (3)  calculates  the  mass  loss  Tate  for  each  wood  component 
(cellulose,  lignin,  and  hemi cellulose)  in  each  slice  using  the  rate  equations, 
(4)  multiplies  these  rates  by  the  local  heat  of  combustion  of  the  volatiles 
generated  by  each  component,  and  (5)  sums  these  contributions  over  the  depth 
of  the  specimen  to  obtain  the  total  heat  release  rate  assuming  complete 
combustion  of  the  volatiles  leaving  the  front  surface.  The  flow  chart  for  the 
computer  program  written  to  Implement  these  calculations  is  shown  in  figure 
2.  The  boundaries  of  each  slice  move  as  the  specimen  shrinks  during  the 
burning  period  so  that  no  solid  material  crosses  a  boundary.  All  of  the 
volatiles  pass  through  the  front  surface.  The  transit  time  of  the  volatiles 
from  their  generation  site  to  the  surface  is  neglected.  However,  the  internal 
surface  convective  heat  transfer  coefficient  is  assumed  to  be  so  large  that 
the  volatiles  are  maintained  in  thermal  equilibrium  with  the  char  through 
which  they  pass.  Secondary  chemical  reactions  with  the  char  layer  are 


neglected.  The  rate  of  ness  generation  for  each  component  la  expressed  as  a 
function  of  the  temperature  and  the  mass  retention  fraction  of  that  component 
for  each  slice  using  an  Arrhenius  type  expression.  The  mass  retention 
fraction  is  the  ratio  of  the  mass  at  any  time  to  the  original  mass.  It  is 
assumed  that  the  Instantaneous  heat  of  combustion  of  the  volatiles  generated 
by.  each  component  depends  only  on  the  mass  retention  fraction  for  that 
component.  This  implies  that  the  chemical  composition  of  each  component  is  a 
unique  function  of  its  mass.  The  temperature  profile  as  a  function  of  time  is 
calculated  from  the  energy  equation  using  finite  differences.  The  rate  of 
change  of  enthalpy,  heat  conduction,  internal  heat  generation  or  absorption 
and  convective  cooling  by  the  flow  of  volatiles  is  taken  into  account.  The 
thermal  properties  are  assumed  to  be  a  function  of  the  temperature  and  the 
total  mass  retention  fraction  for  each  slice.  Then  the  Instantaneous  mass 
retention  fraction  profile  is  determined  by  accounting  for  the  cumulative  loss 
of  volatiles.  An  adiabatic  boundary  condition  is  assumed  at  the  rear  surface 
while  the  front  face  is  exposed  to  a  constant  external  radiant  flux.  Both  of 
these  conditions  can  be  made  more  general  later  if  desired.  The  front  surface 
also  loses  heat  by  reradiatlon.  Up  to  the  time  of  ignition  it  also  loses  heat 
by  laminar  free  convection.  This  is  Justified  for  small  surfaces  such  as 
those  in  the  heat  release  rate  calorimeters.  For  wall  fires  it  will  be  neces¬ 
sary  to  take  turbulence  into  account.  After  ignition  the  convective  heat 
transfer  from  the  flame  to  the  surface  is  taken  into  account  along  with  flame 
radiation.  The  convective  heat  transfer  coefficient  is  multiplied  by  a 
blowing  factor  to  account  for  the  effect  of  the  volatiles  leaving  the  surface. 
Ignition  is  assumed  to  occur  when  the  calculated  total  heat  release  rate 
reaches  30  kW/m^  which  is  the  minimum  heat  release  rate  required  to  maintain  a 
flaw  on  the  surface.  The  use  of  this  ignition  criterion  will  be  discussed  in 


I 


section  7.  The  flame  is  extinguished  when  the  calculated  heat  release  rate 
again  drops  below  30  kW/m^  at  the  end  of  the  flaming  phase.  Only  the  heat 
release  rate  during  the  flaming  phase  is  treated  by  this  model. 


The  mass  retention  fraction  Z^  for  slice  i  is  defined  by 


V"?  ■  f-  lx  ‘y  lz 
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(2-1) 


where  m^  is  the  mass  of  slice  1  on  an  oven  dry  basis,  p  is  its  density,  and 

1,1,  and  l  are  the  contraction  factors  (l.e.  the  ratios  of  the  thickness, 
x  y  z 

width  and  height  of  the  slice  to  their  original  values).  The  superscript  and 
subscript  o  refer  to  the  initial  values. 


The  contributions  of  the  individual  wood  cong>onents  (e. g.  cellulose, 
hemicellulo8e  and  lignin)  are  additive  so  that 


(2-2) 


where  Zj^  is  the  mass  of  the  k  th  component  of  slice  1  divided  by  the 
original  mass  of  the  whole  slice.  It  is  necessary  to  know  the  reaction  rate 
and  heat  of  combustion  of  the  volatiles  for  each  component  as  a  function  of 
the  temperature,  T^,  and  the  component  mass  retention  fraction,  ^/Z®. 

Here  z£  is  original  mass  fraction  of  the  k  th  component. 


The  moisture  is  taken  into  account  by  assigning  a  moisture  retention 
fraction 


WW 


where  hl.  Is  Che  mass  of  adsorbed  wacer  located  in  slice  1  at  any  cioe.  The 
1 

water  produced  by  cheaical  decomposition  Is  accounted  for  In  the  Z  term. 

The  mass  loss  rate  on  a  dry  basis  and  heat  release  rate  per  unit  area  of 
the  original  specimen  are  calculated  by  the  following  formulas: 
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where  Zz  s  Z(t)  and  AXq  Is  the  original  thickness  of  each  slice  (assumed  to  be 
Identical) 


C2-5) 


where  the  heat  of  combustion  of  the  volatiles  Is  given  by 


fJ  -  r  zj 


(2-6) 


since  -  Z±(t)  and  Z^-  Zi>k(t)  then  ■  Ft(t). 


The  energy  equation  Is  first  solved  for  the  temperature.  Then  the  rate 
equations  are  used  to  update  the  mass  retention  fraction  for  each  component  In 
each  slice.  The  total  mass  retention  fraction,  Z^,  for  each  slice  is  then 
determined,  the  property  values  are  updated  and  the  energy  equation  is  solved 
for  the  tetq>erature  at  the  next  time  step. 
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The  Increase  In  enthalpy  of  an  interior  slice  is  equal  to  the  heat 
conducted  in  tainus  the  heat  conducted  out  alnus  the  heat  absorbed  by  the 
pyrolysis  gases  passing  through  on  their  way  to  the  front  surface  minus  or 
plus  the  heat  absorbed  or  generated  by  pyrolysis  and  evaporation.  The  energy 
equation  for  an  interior  slice  can  be  written  as 
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Assuming  instantaneous  mass  transfer. 
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Here  fc"  refers  to  the  mass  flow  per  unit  original  area  of  the  specimen,  not 
the  contracted  area.  The  heat  of  pyrolysis  of  the  slice,  hp^,  is  equal  to  the 
weighted  averages  of  the  heat  of  pyrolysis  of  each  component. 


Since  Zt+At  is  not  known  ahead  of  the  calculation  of  the  temperature  the 
time  step  must  be  small  enough  to  ensure  that  the  approximations 
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are  valid.  In  the  storage  term  (l.h.s.  of  equation  (2-7))  it  is  easier  to 
consider  the  heat-up  of  all  of  the  mass  present  in  the  slice  at  the  beginning 
of  the  time  step  including  that  lost  by  pyrolysis  and  evaporation  during  the 
time  step.  This  mass  loss  then  must  not  be  added  to  the  pyrolysis  gases 
flowing  through  the  slice  as  far  as  transpiration  cooling  of  that  layer  is 
concerned.  Thus  the  mass  flow  rate  through  each  slice  is  uniform,  and  equal 
to  4J+1. 


When  the  energy  equation  is  modified  to  take  the  two  approximations  in 
equations  (2-11)  and  (2-12)  into  account  and  is  divided  through  by  pQ  AXqC^  it 
can  be  rewritten  as 
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It  is  noted  that  AX  -  t  AX  • 
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and  rewrite  the  energy  equation  in  matrix  format  for  solution  of  a  set  of 
algebraic  equations  for  T^+At  using  an  implicit  method: 
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There  is  no  mass  flow  of  volatiles  through  the  rear  face  which  forms  the 


rear  boundary  of  slice  N. 


-  phtS-“*  (p»  *  zl 


Hence  equation  (2-15)  for  i«N  reduces  to 
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Radiative  and  convective  heat  transfer  must  be  taken  into  account  at  the 
front  surface  which  forms  the  front  boundary  of  slice  1.  Absorption  and 
emission  of  radiation  within  the  depth  of  the  specimen  is  assumed  to  be  negli 
gible.  Prior  to  ignition  the  surface  is  cooled  by  laminar  free  convection. 
After  ignition  it  is  heated  by  the  flame.  The  front  surface  temperature,  Ts, 
for  these  calculations  is  obtained  by  linear  extrapolation  from  the  center 
points  of  the  two  slices  nearest  the  surface. 
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The  reradlation  from  the  surface  is  given  by 
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The  energy  flow  across  the  front  surface  is  given  by 
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and 


B  -  m"  (h  +  hf )/Cg 


where  h  Is  Che  convective  heat  transfer  coefficient  from  the  surface  to  the 
ambient  air  which  drops  to  zero  after  ignition,  hf  is  the  convective  heat 
transfer  coefficient  for  the  flame  to  the  surface  which  is  equal  to  zero  prior 
to  ignition,  is  the  flame  temperature,  m"  is  the  mass  flux  through  the 
front  face,  Cg  is  the  heat  capacity  of  the  gas,  and  B  is  the  mass  transfer 
number  introduced  by  Spalding  [14],  which  is  proportional  to  the  mass  loss 
rate  and  thus  accounts  for  blowing.  Prior  to  ignition  the  radiant  heat  flux 
from  the  flame,  q^R,  is  zero. 


The  energy  equation  for  the  first  slice  is  given  by  a  slight  extention  of 
equation  (2-7) 
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In  final  form  Chls  becomes 
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The  change  in  mass  retention  fraction  for  each  component  of  slice  i  during  the 
time  step  is  given  by 
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where  Z^  ^  is  Che  residual  mass  fraction  of  component  k  when  pyrolysis  is 
completed  in  an  inert  atmosphere,  and  and  E^  are  the  frequency  factor  and 
activation  energy  for  component  k.  The  total  mass  retention  fraction  for 
slice  i  is  then  given  by 
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The  mass  retention  fraction  for  water  is  given  by 

-  *  (1  -  Aw  exp  (-  Ew/RTj+At)  At)  (2-25) 

The  parameters  —  K»  p,  C,  C^f  Cg,  a,  Ajj,  A^,  E^,  E^,  F^,  hp^.  ^*v*  ^x*  ^y* 

1  -  are  needed  as  inputs  to  the  model*  While  some  of  these  data  can  be 

z 

obtained  within  the  accuracy  needed  by  consulting  the  literature,  the  thermo- 
physical  properties  K,  p,  C,  a,  i^,  i^,  and  and  the  thermochemical 


material  of  concern.  Most  of  these  parameters  vary  with  temperature  and  the 
mass  retention  fraction. 


3.  THERMOCHEMICAL  PROPERTIES 


3.1  Pyrolyser 


The  pyrolyser  shown  in  figure  3  was  constructed  to  obtain  the  kinetic 


constants  A^  and  E^  and  the  heat  of  combustion  required  by  the  model.  A 
small  cavity  is  created  in  the  center  of  the  aluminum  block  by  bolting  three 
sections  together.  The  block  is  heated  to  a  constant  and  uniform  temperature 


by  a  hot  plate  and  Is  Insulated  by  50  om  thick,  fiberglass  batts.  After  a 
constant  temperature  is  reached,  a  SO  atm  diameter  1  mm  thick  specimen  is 
inserted  by  placing  it  in  the  circular  cutout  in  the  sliding  aluminum  plate 
shown  in  figure  4.  When  the  plate  is  slid  into  position  it  effectively  seals 
off  the  insertion  channel  except  for  a  small  groove  through  which  a  10  mil 
chromel  alumel  thermocouple  passes.  There  are  also  notches  in  the  circular 
portion  of  the  plate  to  allow  easy  flow  of  nitrogen  into  the  specimen  chamber. 
The  thermocouple  is  attached  to  the  specimen  by  threading  it  through  a  series 
of  holes  and  pulling  the  junction  tight  against  the  surface.  A  nitrogen 
stream  sweeps  the  volatile  pyrolysis  products  up  into  a  lean  methane/oxygen 
flame*  The  Increase  in  CO 2  and  H2O  and  the  drop  in  O2  measured  by  the 
analyzers  are  used  to  calculate  the  carbon,  hydrogen,  and  oxygen  content  of 
the  volatiles  according  to  the  procedures  outlined  in  Appendix  C.  The  mass 
loss  rate  necessary  to  calculate  the  kinetic  constants  Ak  and  E^  is  then  given 
by 

fe  -  fcc+  9^+  h0  (3-1) 

The  heat  release  rate  is  calculated  by  determining  the  external  oxygen 
requirement  for  complete  combustion  to  water  and  carbon  dioxide  and  multiply¬ 
ing  by  13100  kJ/kg.  It  takes  8  grams  of  oxygen  to  convert  1  gram  of  hydrogen 
to  water  and  8/3  grams  of  oxygen  to  convert  1  gram  of  carbon  to  carbon 
dioxide.  The  mass  of  oxygen  present  in  the  volatiles  oust  be  subtracted  from 
the  total  required  to  completely  oxidize  the  hydrogen  and  carbon  in  order  to 
determine  the  external  oxygen  requirement.  It  is  assumed  that  13100  kJ  of 
heat  is  released  for  each  kilogram  of  oxygen  consumed  [15].  This  approxima¬ 
tion  is  good  to  within  ±  52  for  most  organic  compounds.  Thus 
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The  gross  heat  of  combustion  of  charred  and  uncharred  specimens  of 
Douglas  fir  particle  board  were  determined  with  the  oxygen  bomb  calorimeter. 
These  data  are  plotted  as  open  circles  in  figure  5.  The  data  point  on  the 
vertical  axis  is  pure  carbon  which  might  be  considered  to  be  the  ideal  end 
point  in  the  charring  process.  However,  in  a  fire  situation  it  is  the  net 
heat  of  combustion,  which  refers  to  the  case  where  the  water  remains  as  a 
vapor  in  the  final  state,  that  is  important.  From  the  overall  chemical  compo¬ 
sition  of  the  wood  the  amount  of  water  produced  during  complete  combustion  can 
be  determined.  Thus  the  heat  of  evaporation  of  the  water  can  be  taken  into 
account  in  order  to  obtain  the  net  heat  of  combustion  of  the  virgin  wood.  The 
upper  dashed  curve  is  the  net  heat  of  combustion  assuming  a  linear  drop  in  the 
correction  factor  for  water  as  the  final  char  is  approached. 


The  heat  of  combustion  of  the  volatiles  Is  equal  to  the  total  heat 
released  during  combustion  of  the  wood  minus  the  total  heat  released  during 
combustion  of  the  char  divided  by  the  difference  in  mass  between  the  wood  and 
the  char.  The  Instantaneous  heat  of  combustion  of  the  volatiles  is  given  by 

P  -  AOftSi  -  AH  +  )  ■  d.ftS.  .  (3-4) 
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The  derived  values  appear  as  closed  circles  in  figure  5. 

Using  a  linear  approximation  to  the  points  between  Z  ■  0. 4  and  1.0  the 
heat  of  combustion  of  the  volatiles  can  be  represented  over  the  range  of 
concern  by 

F  -  20  [1.24  -  Z]  MJ/kg  Z  >  0.26  (3-5) 

-  33.6  Z  <  0.26 

which  is  expressed  by  the  lower  dashed  curve.  Although  the  calculated  points 
demonstrate  some  curvature,  the  uncertainty  in  the  assumed  correction  factor 
for  water  does  not  justify  the  inclusion  of  additional  terms  in  eq.  (3-5). 

3.3  Calculation  of  Kinetic  Constants  from  Mass 
Loss  Data  at  Constant  Temperature 

Although  there  are  several  wood  components  present,  each  with  different 
kinetic  constants,  it  was  decided  to  try  to  fit  the  decomposition  of  particle 
board  to  a  single  first-order  reaction.  This  is  the  normal  assumption  for  the 
mass  loss  rate  models.  In  that  case  for  specimens  at  a  constant  temperature 
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so  that 


(a  -  a  ) 

ln  fao’-"af')  *  *  c  A  «*P  (“  E/RT) 


and 

E  ! 

In  A  -  -  i  •  (3-7). 

The  final  char  fraction  Of/niQ  was  found  to  be  0.225  after  heating  a 
specimen  wrapped  ln  aluminum  foil  ln  an  oven  for  5  minutes  at  560°C.  The 
intercept  and  slope  of  the  plot  of  the  l.h.s.of  equation  (3-5)  against  1/T 
yields  a  frequency  factor  A  of  5.9  x  107  s”1  and  an  activation  energy  Ea  of 
121  kJ/mol  for  the  Douglas  fir  particle  board,  as  seen  in  figure  6.  The 
constancy  of  the  l.h.s.  of  equation  (3-7)  for  a  given  value  of  T  depends  on 
the  assumption  of  a  single  first  order  reaction.  Since  this  assumption  is  not 
completely  valid,  the  expression  in  the  brackets  will  vary  for  different 
combinations  of  m  and  t  even  at  a  fixed  temperature  and  thus  could  be  respon¬ 
sible  for  some  of  the  scatter  in  figure  6.  While  the  data  need  to  be  analyzed 
ln  terms  of  more  than  one  reaction  the  constants  determined  in  figure  6  are 
used  ln  the  initial  checkout  of  the  computer  model  for  heat  release  rate. 


3. 4  Contraction  Factors 


Tha  contraction  factors  measured  along  the  surface  (i.e. ,  it  and  t  )  for 

y  z 

southern  pine  are  plotted  in  figure  7.  The  contraction  is  seen  to  be  ouch 
greater  across  the  grain  than  it  la  along  the  grain.  In  addition,  the  thick¬ 
ness  of  each  specimen  was  measured  with  a  micrometer.  There  was  considerably 
more  scatter  in  the  latter  data  but  they  tended  to  center  around  the  lower 
curve.  The  contraction  factors  for  Douglas  fir  particle  board  are  displayed 
in  figure  8.  The  contraction  is  Independent  of  direction  and  is  midway 
between  the  contraction  factors  for  southern  pine.  The  data  for  Douglas  fir 
particle  board  can  be  represented  by  the  expressions, 

4  -  4  -  4  -  1  for  0.65  <  Z  <  1.0  (3-8) 

x  y  z  —  — 

l  -  1  -  4  -  1  -  (0.65  -  Z)2  for  Z  <  0.65  (3-9) 

x  y  z 

3.5  Heat  of  Pyrolysis 

Vfhlle  differential  scanning  calorimeter  (DSC)  data  has  been  taken  on  both 
southern  pine  and  particle  board  to  determine  the  heat  of  pyrolysis,  these 
data  have  not  been  fully  analyzed  yet.  The  heat  of  pyrolysis  is  suspected  to 
be  small  and  has  been  therefore  set  equal  to  zero  in  the  model.  Atreya  [13] 
has  found  that  the  results  of  his  mass  loss  model  for  wood  were  relatively 
insensitive  to  the  heat  of  pyrolysis.  That  sensitivity  was  also  checked  on 
this  project  by  making  a  run  with  the  model  assuming  a  heat  of  pyrolysis  of 
400  kJ/kg  which  is  approximately  the  value  for  alpha  cellulose.  The  peak  heat 
release  rate  was  reduced  by  20Z. 


3.6  Parameters  Dealing  with  Moisture 


Atreya  found  his  mass  loss  rate  model  to  be  very  sensitive  to  the 
moisture  content  [13].  He  determined  the  frequency  factor  and  activation 
energy  for  the  desorption  of  water  experimentally  to  be  4.5  x  10J  s“*  and 
43.9  kJ/mole.  The  heat  of  vaporization  of  water  is  2.4  MJ/kg.  These  values 
were  used  in  the  present  model. 

4.  THERMOPHYSICAL  PROPERTIES 

4. 1  Heat  Capacity 

The  specific  heat  of  wood  does  not  depend  on  density  or  even 
significantly  on  the  species.  However,  it  does  depend  rather  strongly  on 
temperature.  In  1912  Dunlap  [16]  described  this  temperature  dependence  up  to 
106°C  for  20  different  species  of  wood  as 

C  -  1.11  +  0.00486  T  (4.1) 

In  1969  Koch  [17]  found  that 

C  -  1.11  +  .0042  T  (4-2) 

for  spruce  pine  between  60°C  and  140°C  where  T  Is  in  °C  and  C  is  in  kJ/kg°K. 
The  present  model  simply  extends  the  linear  relationship  In  equation  (4-2)  to 
the  highest  temperature  at  which  the  wood  retains  its  chemical  structure.  Fo 
moist  specimens  the  heat  capacity  of  the  wood  and  the  heat  capacity  of  the 


moisture  are  linearly  combined.  The  heat  capacity  of  charcoal  is  reported  to 
be  0.67  kJ/kg  «K  [16].  For  the  purposes  of  the  initial  running  of  the  model 
the  heat  capacity  has  been  assumed  to  vary  linearly  between  wood  and  final 
char  and  to  have  the  same  temperature  dependence  as  for  wood.  Thus 

C  -  (1.11  +  0.0042  T)  (0.43  +  0.57  Z)  (4-3) 

for  oven-dry  wood  and  char. 

4.2  Thermal  Conductivity  Data  from  the  Literature 

A  universal  formula  for  the  thermal  conductivity  of  wood  (normal  to  the 
grain)  and  plywood  (across  the  thickness)  as  a  function  of  density  and  mois¬ 
ture  content  was  established  experimentally  by  MacLean  [18].  This  is  given  by 

K  -  [237  +  2.00  p  (1  4  .020  M)]  x  10-7  kW/m.K.  (4-4) 

The  density,  p,  is  expressed  in  kg/m^  over  a  range  of  300-800  kg/nr*  and  the 
moisture  content,  m,  in  percent  by  weight  over  a  range  of  0-40Z.  This  formula 
was  also  found  to  hold  for  charred  wood  up  to  a  37Z  mass  loss.  It  is  noted 
that  when  p  goes  to  zero,  k  is  equal  to  the  thermal  conductivity  of  air. 

The  ratio  of  the  thermal  conductivity  along  the  grain  to  that  across  the 
grain  will  depend  on  the  fibril  angle.  A  typical  ratio  is  2.5.  The  more 
nearly  aligned  the  fibers  are  along  the  length  of  the  cell  wall  the  higher  the 
thermal  conductivity  in  grain  direction. 


Thermal  conductivity  data  on  oven  dry  flberboard  and  particle  board  at 
24°C  was  obtained  by  Lewis  [19]  and  is  presented  in  Table  1. 

Using  these  data  and  assuming  that  the  contributions  of  the  water  and  air 
were  the  same  in  the  flberboard  and  particle  board  as  they  were  in  the  wood, 
equations  of  the  same  form  as  equation  (4—3 )  were  obtained  for  these 
materials.  For  flberboard 

K  -  [237  +  1.11  p  (1  *  .036  M)]  x  10"7lcW/m.K  (4-5) 

and  for  particle  board 

K  -  [237  +  1.37  p  (1  4  .029  M)]  x  10"7kW/m.K.  (4-6) 

Lewis  also  reported  average  Increases  in  thermal  conductivity  of  1.37  x  10~7 
and  2.41  x  10“7  per  deg  C,  respectively,  for  flberboard  and  particle  board. 

The  temperature  dependence  of  the  thermal  conductivity  of  wood  was  given  by 
Kollmann  [20]  as 

K  -  Kq  [1  +  (1.1  -  9.8  x  10"%)  (T  -  Tq)/100]  kW/m.K  (4-7) 

The  upper  limit  for  this  formula  is  100 °C. 


4.3  Thermal  Conductivity  Measurements 


The  thermal  conductivity  was  measured  on  this  project  as  a  function  of 
temperature  up  to  2S0°C  by  exposing  vertical  specimens  of  thicknesses  between 
1  and  2  mm,  which  are  at  least  40  mm  wide  and  high,  to  a  radiant  panel. 
Thermocouples  were  attached  to  the  front  and  rear  surfaces.  At  equilibrium 
the  heat  flowing  through  the  specimen  is  lost  from  the  rear  face.  Thus 


K  (Tj  -  T2)/D  -  ao  (T*  -  T*)'  +  h2  (T2  -  Tq) 


[«  (T2  "  To)  *  h2  (h  -  *.)] 


D/CTj  -  T2) 


(4-8) 


where  T ^  and  T2  are  the  front  and  rear  surface  temperatures,  respectively. 

The  only  unknown  on  the  r.h.s.  (i.e. ,  t^)  could  be  calculated  assuming  laminar 
free  convection.  Thus  according  to  Schlichting  [21], 


h  -  0.636  ,  “/*  (^)1M 


(4-9) 


where  x  is  the  distance  from  the  bottom  of  the  specimen  to  the  thermocouple. 


It  can  also  be  calculated  from  the  energy  balance: 


*1  ■  ««  (Tt  *  A.  *  2To)  *  hl  (T1  *  To)  +  h2  (T2 


(4-10) 


If  the  temperature  difference,  T2  -  Ti,  is  small  enough  that  X  and  v  are  the 
same  at  both  faces  and  the  thermocouple  junctions  are  at  the  same  height,  the 
ratio  hj/h2  is  found  from  equation(4-9)  to  be 


As  seen  in  figure  9  the  5  mil  (0.127  mm)  diameter  thermocouples  attached 
to  the  front  and  rear  surface  provide  support  for  the  specimen.  The  junctions 
are  held  firmly  against  the  surface  by  threading  them  through  small  diameter 
holes  and  pulling  ^lght.  Improved  thermal  contact  as  well  as  a  high  and  known 
absorptivity  is  achieved  by  coating  both  surfaces  of  the  specimen  with  a  spray 
coat  of  Krylon  1602*  which  has  an  absorptivity  of  0.94  [22].  The  effect  of 
the  extra  thickness  provided  by  the  paint  has  not  been  investigated  yet.  The 
wires  then  go  to  a  digital  thermocouple  readout  with  a  built-in  electronic 
cold  junction. 

The  external  radiant  flux  incident  on  the  front  face  is  controlled  by 
varying  the  distance  from  the  panel.  Equilibrium  is  established  in  less  than 
one  minute  thereby  allowing  the  specimen  to  go  to  relatively  high  temperatures 
without  significant  decomposition  as  opposed  to  the  standard  guarded  hot  plate 
method  which  takes  many  hours  to  establish  thermal  equilibrium.  At  each 
distance  the  external  heat  flux  and  surface  temperatures  are  recorded.  A 

*The  product  name  is  given  for  identification  purposes  only  and  does  not 
constitute  an  endorsement. 


snail  conpuCer  program  was  written  to  calculate  the  thermal  conductivity  from 
the  above  input  data.  The  convective  heat  transfer  coefficient  is  also  calcu¬ 
lated  by  the  program  as  a  check  on  the  system. 

Table  2  shows  the  results  for  a  1.27  mm  thick  oven  dry  specimen  of 
southern  pine  while  table  3  shows  similar  data  for  a  charred  specimen.  These 
data  are  plotted  in  figure  10.  The  lines  represent  a  least  squares  fit.  Ihe 
extrapolated  value  of  the  thermal  conductivity  at  20°C  is  1.22  x  10-4 
kW/m  »K  compared  to  1.23  x  10-4  kW/m  »K  predicted  by  MacLean's  formula  (4-4) 
using  the  measured  density  of  southern  pine  which  is  498  kg/nr*.  The 
calculated  values  from  eq.  (4-4)  for  both  charred  and  uncharred  pine  are 
included  in  the  figure. 

The  measured  density  of  the  char  is  207  kg/m*.  However,  the  calculated 
temperature  coefficient  using  Kollman's  formula  (4-7)  is  an  order  of  magnitude 
higher  than  the  measured  one.  The  thermal  conductivity  data  plotted  in  figure 
10  can  be  expressed  by  the  formula 

K  -  1.20  (1  +  0.0026  (1.35  -  Z)  T)  (0.71  Z  +  0.29)  x  10"4  ktf/m.K  (4-13) 
where  the  temperature  is  in  °C. 

The  experimentally  determined  thermal  conductivity  of  uncharred  and 
charred  Douglas  fir  particle  board  is  listed  in  tables  4  and  5  and  plotted  in 
figure  11.  The  value  of  the  thermal  conductivity  of  uncharred  particle  board 
at  20°C  in  figure  11  was  determined  with  a  guarded  hot  plate.  The  upper  line 
was  obtained  by  a  least  squares  fit.  Because  of  the  larger  scatter  and  the 
limited  temperature  range  of  the  data  the  use  of  a  least  squares  fit  for  the 


lower  line  led  to  an  unreasonably  high  value  of  the  thermal  conductivity  of 
the  char  at  ambient  temperature.  Thus  the  lower  line  was  passed  through  the 
center  of  gravity  of  the  data  points  and  the  calculated  value  of  the  thermal 
conductivity  using  eq.  (4-6)  with  the  measured  char  density.  The  thermal 
conductivity  represented  by  the  lines  in  figure  11  can  be  expressed  as 

K  -  1.16  (1  +  0.0058  (1.6  -  Z)  T)  (0.84  Z  +  0.16)  x  10_4  kW/m.K  (4-14) 

Equation  (4-14)  was  used  in  the  computer  model.  This  equation  was  determined 
for  Che  temperature  range  over  which  wood  could  be  expected  to  exist  in  its 
original  form  and  is  simply  extrapolated  into  the  high  temperature  region. 
High  temperature  data  on  charred  specimens  are  needed  to  accurately  predict 
the  thermal  conductivity  over  the  full  temperature  range  required.  The 
measured  value  of  the  thermal  conductivity  at  20°C  using  eq.  (4-14)  for 
extrapolation  is  1.24  x  10“^  kW/mK  compared  to  1.21  x  10“^  kW/m  »K  calculated 
from  equation  (4-6)  using  a  density  of  709  kg/m^  for  Douglas  fir  particle 
board.  Ihe  calculated  and  measured  temperature  coefficients  were  within  20% 
for  the  particle  board. 

The  rear  surface  convective  heat  transfer  coefficient  was  calculated  for 
comparison  with  the  measured  values  reported  in  table  2.  At  the  highest  heat 
flux  the  rear  surface  temperature  for  uncharred  southern  pine  was  234°C.  In 
this  case  the  average  film  temperature  in  the  boundary  layer  was  127 ®C.  The 
thermal  conductivity  and  kinematic  viscosity  at  this  temperature  are 
3.37  x  10”5  kW/mK  and  25.9  x  10~^m^/s  [21],  respectively.  The  thermocouple 
junction  was  0.02  m  above  the  lower  edge  and  the  ambient  temperature  was 
30*C.  When  these  data  were  substituted  into  equation  (4-9)  the  heat  transfer 


than  the  measured  value.  In  addition  to  experimental  errors,  the  flow  may  not 
correspond  to  the  ideal  condition  required  by  equation  (4-8).  On  the  other 
hand  the  height  of  the  particle  board  specimens  was  twice  that  of  the  pine 
yielding  a  calculated  coefficient  at  230°C  of  0*0079  which  is  30%  higher  than 
the  measured  value  in  table  4.  These  30%  differences  only  lead  to  a  10%  error 
in  the  thermal  conductivity  at  230 °C. 

5.  CALCULATIONS  OF  HEAT  RELEASE  SATE  USING  THE  MODEL 

The  computer  program  requires  input  data  for  a  number  of  parameters. 

Some  of  these  have  been  experimentally  determined  on  the  project,  some  have 
been  taken  from  the  literature  and  some  are  simply  assumed  at  this  stage. 

These  parameters  along  with  their  assigned  values  and  sources  are  listed  in 
Table  6.  The  front  and  rear  surface  temperature,  the  mass  loss  rate  (equation 
(2-4)),  the  heat  release  rate  (equation  (2-5)),  the  char  depth,  and  the  effec¬ 
tive  heat  of  combustion  are  all  calculated  as  a  function  of  time  during  the 
test.  The  effective  heat  of  combustion  is  simply  equal  to  the  ratio  of  the 
heat  release  rate  and  the  mass  loss  rate.  The  char  depth  is  equal  to  AXq 
multiplied  by  the  number  of  slices  for  which  Z  is  less  than  0.95.  Some  of 
these  outputs  are  plotted  in  figures  12  to  19  for  different  external  radiant 
heat  fluxes,  specimen  thicknesses,  and  moisture  contents.  The  baseline  case 
is  for  an  oven  dry  12.7  mm  thick  specimen  of  Douglas  fir  particle  board 
exposed  to  an  external  radiant  heat  flux  of  50  kW/m^.  In  all  cases  the  rear 
surface  is  assumed  to  be  insulated  and  in^ervlous  to  the  flow  of  volatile 
pyrolysis  products  and  water  vapor. 


The  program  has  been  set  up  on  a  Tektronix  4052  terminal  using  BASIC. 

The  total  number  of  slices  is  limited  to  60  by  the  capacity  of  the  micro¬ 
computer.  The  longest  acceptable  time  step  is  one  second  at  an  incident  flux 
of  25  kW/m^,  and  0. 1  second  at  an  incident  flux  of  100  kW/m^.  This  can  be 
increased  later  on  in  the  run.  These  critical  times  were  determined  by  using 
successively  smaller  time  steps  until  there  were  no  further  changes  in  the 
output  data  curves.  The  calculations  require  about  30  seconds  per  time  step. 

A  typical  run  takes  about  4  hours. 

The  calculated  front  and  rear  surface  temperatures  for  a  dry  12.7  mm 
specimen  of  Douglas  fir  particle  board  exposed  at  an  external  radiant  flux  of 
50  kW/er1  are  shown  in  figure  12.  The  front  surface  temperature  shows  a  cusp 
at  45  seconds,  when  ignition  occurred.  It  approaches  700°C  by  the  time  the 
flame  goes  out.  Figure  13  shows  the  predicted  progression  of  the  char  front 
into  the  material.  The  char  front  here  is  taken  to  be  located  at  the  center 
point  of  the  slice  whose  mass  retention  fraction  has  just  dropped  to  0.95. 
After  an  Initial  rapid  penetration  the  charring  rate  becomes  nearly  constant 
until  the  rear  surface  temperature  starts  to  rise.  Then  it  Increases 
abruptly.  The  apparently  constant  rate  may  depend  on  specimen  thickness. 
Figure  14  shows  the  predicted  mass  flow  through  the  front  surface.  The  second 
large  peak  occurs  when  the  thermal  wave  is  reflected  from  the  insulated  rear 
face.  This  second  peak  would  be  missing  altogether  if  the  rear  surface  were 
maintained  at  some  low  temperature.  Figure  15  shows  the  calculated  heat 
release  rate  curve.  Its  features  are  similar  to  the  mass  loss  rate  curve 
except  that  it  rises  abruptly  from  zero  to  30  kW/ts^  at  ignition  and  the  ratio 
of  the  second  and  first  peak  is  higher  for  the  heat  release  rate  curve.  The 
first  peak  is  about  240  kW/m^.  The  effective  heat  of  combustion  is  plotted  in 
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figure  16«  It  is  relatively  constant  over  the  first  part  of  the  curve  decay¬ 
ing  only  slightly  to  a  minimum  and  then  rising  abruptly  when  the  rear  surface 
begins  to  char.  The  magnitude  during  the  early  part  of  the  curve  is  12  MJ/kg 
versus  20  MJ/kg  for  the  gross  heat  of  combustion  of  the  uncharred  particle 
board  determined  by  the  oxygen  bomb  calorimeter.  The  bump  in  the  curve  around 
400  s  is  an  artifact  caused  by  a  change  in  the  time  step  size  at  that  point. 

Throughout  the  above  calculations,  the  heat  of  pyrolysis  has  been  set 
equal  to  zero.  The  effect  of  the  heat  of  pyrolysis  on  the  heat  release  rate 
curve  is  displayed  in  figure  17  where  it  has  been  set  equal  to  400  kJ/kg  which 
is  approximately  the  value  for  alpha  cellulose.  This  resulted  in  a  20Z  drop 
in  the  peak  heat  release  rate  and  a  substantially  larger  drop  in  the  second 
peak.  Figure  18  shows  the  effect  of  a  7Z  moisture  content.  This  is  the 
equilibrium  value  for  Douglas  fir  conditioned  at  50Z  RH  at  23 °C.  This  results 
in  a  drop  in  the  peak  heat  release  rate  of  about  10Z.  The  value  of  the  first 
peak  is  relatively  unaffected  by  decreasing  the  thickness  of  the  slab  but  the 
second  peak  is  much  higher  and  earlier  in  time  as  seen  in  figure  19.  Lowering 
the  irradiance  to  25  kW/m^  lowers  the  first  peak  and  extends  the  time  scale  as 
seen  in  figure  20.  The  calculated  heat  release  rate  curves  for  75  and  100  kW 
are  presented  in  figures  21  and  22. 

6.  HEAT  RELEASE  RATE  MEASUREMENTS  IN  THE  CONE  CALORIMETER 

Heat  release  rate  measurements  were  made  on  oven  dry  specimens  of  Douglas 
fir  particle  board  at  external  radiant  flux  levels  of  25,  50,  75,  and 
100  kW/m The  specimens  were  all  12.7  mm  thick.  The  measurements  were  made 
in  the  vertical  orientation  in  the  NBS  cone  calorimeter  [7]  and  are  shown  in 
figures  23  through  26. 


The  hut  release  races  are  determined  using  Che  scandard  method  for  this 


calorimeter  which  is  based  only  on  Che  oxygen  concentration  and  the  volume 
flow  measured  in  the  exhaust  duct.  A  standard  estimated  correction  is  made  to 
account  for  the  difference  between  the  molar  flow  of  air  into  the  system  and 
the  molar  flow  of  gas  in  the  exhaust  duct  due  to  chemical  reaction.  The  mass 
loss  rate  is  determined  from  the  load  cell. 

The  first  peaks  of  the  measured  heat  release  rate  curves  increase  with 
the  Incident  radiant  flux  but  are  about  20Z  lower  than  the  calculated  ones. 

The  second  peaks  are  lower  than  the  first  peaks  for  the  measured  curves.  This 
is  probably  due  to  the  absorption  of  some  heat  by  the  insulation  material  in 
contact  with  the  rear  surface.  The  first  peaks  fall  off  faster  and  the  second 
peaks  occur  later  for  the  calculated  curves.  This  discrepancy  could  be  due  to 
Inadequate  data  on  the  thermal  properties  of  the  char  at  high  temperature  or 
unrealistic  kinetics  defined  by  the  assumption  of  a  single  first  order 
reaction.  The  calculated  heat  of  combustion  is  close  to  the  measured  values 
and  la  Independent  of  the  incident  radiant  flux.  See  for  example  figures  16 
and  27-30.  The  early  segments  of  these  experimental  curves  to  the  left  of  the 
dashed  lines  are  due  to  instrumental  problems  and  should  be  ignored. 

The  three  analyzer  method  of  calculating  the  heat  release  rate  and  effec¬ 
tive  heat  of  combustion  in  the  cone  calorimeter  relies  on  the  concentrations 
measured  in  the  carbon  dioxide  and  water  vapor  analyzers  as  well  as  the  oxygen 
analyzer. 


The  density  of  the  exhaust  gases  is  determined  from  the  concentrations  of 
oxygen,  carbon  dioxide  and  water  vapor  in  the  analyzers  according  to  Appendix 


A.  This  density  along  with  the  exhaust  temperature  and  pressure  drop  across 
the  orifice  plate  in  the  duct  determine  the  nitrogen  flow  through  the  system 
according  to  Appendix  B.  The  nitrogen  flow  along  with  the  above  concentra¬ 
tions  were  used  in  the  calculation  of  the  mass  flows  of  carbon,  hydrogen  and 
oxygen  from  the  specimen  using  Appendix  C.  The  latter  quantities  were  used  to 
determine  the  mass  loss  rate  and  heat  release  rate  and  hence  the  effective 
heat  of  combustion  using  the  methods  of  section  3.2.  The  calculation  of  heat 
release  rate  using  all  three  analyzers  is  based  on  the  assumption  that  the 
fuel  vapors  generated  by  the  burning  specimen  are  composed  only  of  carbon, 
hydrogen,  and  oxygen  and  that  combustion  is  complete.  It  is  not  recommended 
for  the  determination  of  heat  release  rates  in  general.  The  heat  release  rate 
and  the  effective  heat  of  combustion  determined  by  the  three  analyzer  method 
are  shown  in  figures  31  and  32  for  the  50  kU/m^  external  radiant  flux.  These 
curves  depart  from  those  determined  by  the  standard  method  early  in  the 
exposure  because  of  as  yet  unresolved  problems  with  the  time  response  of  the 
water  vapor  measurement  system. 

7.  DISCUSSION 

The  prediction  of  the  heat  release  rate  of  wood  involves  (1)  the  develop¬ 
ment  of  a  suitable  computer  program  which  includes  all  of  the  Important 
physics  with  acceptable  assumptions  and  approximations  and  which  converge  to 
the  exact  solution  of  the  formulated  equations  and  (2)  the  gathering  of  the 
required  input  data  on  the  properties  of  the  particular  wood  specimen  over  a 
large  temperature  and  decomposition  range. 


The  coaputer  model  Includes  a  description  of  the  heat  transfer  at  the 
surfaces,  an  ignition  criterion,  the  foraulation  of  an  energy  equation  in  the 
interior  along  with  rate  equations  for  the  various  chemical  species  which  make 
up  the  solid  wood,  a  means  of  accounting  for  the  shrinkage  of  the  char,  and  a 
way  of  dealing  with  the  fissures  formed  by  the  shrinkage.  This  model  only 
applies  during  the  flaming  phase.  It  is  assumed  that  during  this  period  the 
oxygen  is  completely  excluded  from  the  surface  so  that  there  is  no  oxidation 
of  the  char  and  the  pyrolysis  takes  place  in  an  inert  atmosphere. 

The  heat  transfer  to  the  surface  includes  a  specified  external  radiant 
flux,  radiation  from  the  surface  and  convective  heat  losses  assuming  laminar 
free  convection.  Laminar  flow  seems  reasonable  for  the  100  mm  high  specimens 
in  the  cone  calorimeter  which  is  used  to  check  out  the  model.  Turbulence  will 
have  to  be  Included  in  dealing  with  the  heat  release  rate  from  a  burning 
wall.  After  ignition  the  convective  loss  term  is  replaced  by  flame  radiation 
and  convective  heat  transfer  to  the  surface.  The  radiation  to  the  surface  was 
assumed  to  be  10  kW/m2  based  on  measurements  by  Tewarson  [8].  This  would  only 
be  expected  to  apply  to  small  surfaces.  The  convective  heat  transfer  coeffi¬ 
cient  was  obtained  by  measuring  the  total  flux  to  a  3  on  diameter  water  cooled 
heat  flux  gage  embedded  in  the  surface.  The  10  kW/m2  flame  radiation  was 
subtracted  from  the  30  kW/m2  measured  by  the  gauge.  This  difference  was 
divided  by  the  difference  between  the  flame  temperature,  which  was  measured  by 
a  2  mil  thermocouple  to  be  1200 *C,  and  the  water  temperature  of  20°C.  This 
yielded  a  convective  heat  transfer  coefficient  of  0.017  kW/m  «K.  Any  tempera¬ 
ture  dependence  was  neglected.  However,  the  coefficients  both  before  and 
after  ignition  are  modified  by  a  blowing  factor.  This  depends  on  the  B  number 
which  is  equal  to  the  product  of  the  mass  loss  rate  times  the  heat  capacity  of 


Che  volatiles  divided  by  Che  appropriate  heat  transfer  coefficient.  The  heat 
capacity  of  the  volatiles  was  assumed  to  be  the  same  as  that  of  air.  Emission 
and  absorption  of  radiation  at  the  surface  of  the  specimen  depends  on  the 
surface  emissivity.  This  was  determined  to  be  0.80  from  the  ratio  of  the 
derivatives  of  the  surface  temperature  with  respect  to  the  square  root  of  the 
time  at  the  beginning  of  the  exposure  for  unblackened  and  blackened  specimens 
exposed  to  a  gas-fired  radiant  panel.  Any  change  in  spectral  emissivity  with 
wave  length  was  neglected.  A  state  of  the  art  flame  heat  transfer  model  needs 
to  be  incorporated  into  the  heat  release  rate  model. 

The  piloted  ignition  criterion  adopted  here  was  the  attainment  of  a 
calculated  heat  release  rate  which  would  be  sufficient  to  maintain  a  flame 
over  the  entire  surface  of  the  specimen.  For  lower  heat  release  rates  the 
shortening  of  the  flame  will  lead  to  extinction.  In  the  analysis  of  the 
Steiner  tunnel  [23]  it  was  found  that  this  critical  value  was  approximately 
50  kW/m2.  However,  prior  to  ignition  when  the  test  calculation  is  made  there 
is  no  flame  heat  transfer.  Hence  the  calculated  heat  release  rate  would  be 
less  than  50  kW/m  at  the  time  just  prior  to  piloted  ignition.  It  was  reduced 
to  30  kW/m2  so  that  ignition  could  be  achieved  in  the  model  at  an  external 
radiant  heat  flux  of  25  kW/m2  which  is  sufficient  to  cause  ignition  of  this 
material  in  the  cone  calorimeter.  This  choice  could  have  considerable  impact 
on  the  time  to  ignition  but  probably  a  minor  effect  on  the  heat  release  rate 
beyond  the  very  early  part  of  the  curve.  In  any  event  this  is  another  part  of 
the  model  which  needs  to  be  improved. 

It  was  necessary  to  use  a  difference  formulation  for  the  energy  equation 
because  the  dimensions  of  the  slices  and  the  material  properties  change  during 


the  burning  period  and  are.  In  general,  different  for  each  slice.  A  grid  of 
60  slices  was  necessary  to  prevent  oscillations  In  the  calculated  heat  release 
rate  versus  time  curve  for  the  12.7  mm  Douglas  fir  specimens.  This  grid  size 
required  a  time  step  no  greater  than  1  second  at  an  external  radiant  flux  of 
25  kW/m^  and  0. 1  second  for  100  ktf/m^.  However,  after  the  minimum  in  the  heat 
release  rate  curve  the  time  step  could  be  increased.  The  maxima  permissible 
time  step  was  found  by  reducing  it  until  there  was  no  perceptible  difference 
in  the  generated  heat  release  rate  curve.  The  time  step  determined  this  way 
did  not  satisfy  the  usual  convergence  criterion, 

-■SA*--,  <  1/2,  (7-1) 
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in  the  case  of  the  larger  time  steps.  If  the  oscillations  referred  to  above 
were  due  to  the  violation  of  the  convergence  criterion,  the  effect  of  reducing 
the  grid  spacing  would  be  to  aggravate  the  problem.  They  are  more  likely  to 
be  due  to  intolerably  large  changes  in  z  or  in  the  material  properties  over 
the  distance  Ax. 

An  important  element  of  this  model  is  the  ability  to  consider  the  rate 
equation,  heat  of  pyrolysis,  heat  of  combustion  of  the  volatiles  generated 
during  pyrolysis,  and  heat  capacity  separately  for  each  chemical  component 
present  in  the  solid  wood.  However,  as  a  practical  consideration  it  is  neces¬ 
sary  to  lump  in  so  far  as  possible  all  of  the  species  with  similar  properties 
and  ignore  those  that  are  present  in  very  small  quantities.  Thus  wood  might 
be  considered  to  consist  of  three  components  -  cellulose,  lignin  and  hemi- 
cellulose  -  although  hemicellulose  consists  of  a  number  of  different  carbo¬ 
hydrates  and  there  is  a  small  quantity  of  miscellaneous  extractives  present. 
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To  make  use  of  the  additional  detail  it  Is  necessary  to  measure  the  properties 
of  the  individual  components.  Since  this  has  not  been  done  as  yet,  it  has 
been  necessary  to  consider  wood  as  a  single  chemical  species  in  this  report. 
Assigning  a  single  first  order  reaction  to  wood  has  not  been  very  satisfactory 
as  can  be  seen  in  figure  6  where  different  combinations  of  m  and  t  at  a  fixed 
temperature  do  not  lead  to  the  same  value  for  the  ordinate  as  they  would  if 
the  rate  were  given  by  equation  (3-6). 

The  char  shrinkage  normal  to  the  surface  is  easily  handled  by  the  energy 
equation  by  taking  the  change  in  slice  dimensions  into  account.  The  reduction 
In  spacing  results  in  higher  temperature  gradients  which  Increase  heat  conduc¬ 
tion.  Shrinkage  parallel  to  the  surface,  while  the  specimen  width  is  main¬ 
tained  by  the  rigidity  of  the  remaining  virgin  wood,  sets  up  stresses  which 
relieve  themselves  by  the  formation  of  fissures.  These  are  visible  as  deep 
cracks  in  the  surface.  The  reduction  in  -the  area  of  the  region  of  the  surface 
between  the  fissures  represents  a  reduction  in  the  area  through  which  heat  Is 
conducted  and  also  In  the  area  normal  to  the  external  heat  flux  and  the  heat 
flux  from  the  flame.  This  results  in  a  reduction  in  heat  transfer.  However, 
the  radiant  flux  can  now  penetrate  deep  within  the  specimen  in  the  area  of  the 
fissures.  This  Is  counteracted  by  the  easy  outflow  of  volatiles  through  the 
fibers  parallel  to  the  surface.  They  flow  into  the  fissures  and  carry  the 
excess  heat  with  them  as  they  exhaust  through  the  cracks.  Up  until  the  time 
of  the  fissure  formation  the  volatiles  cool  the  char  layer  as  they  pass 
through  it.  Then  the  cooling  is  shut  off  as  the  layer  is  by-passed.  The 
depth  of  the  fissures  is  taken  to  be  the  middle  of  the  slice  whose  mass  reten¬ 
tion  factor  has  just  dropped  to  Zf  +  0.05  (1  -  Zf)»  This  is  where  95Z  of  the 


mass  that  can  be  lost  is  lost 


The  bulk  properties  required  as  Input  to  the  model  Include  the  thermal 
conductivity,  heat  capacity,  char  contraction  factors,  kinetic  constants,  and 
the  heat  of  combustion  and  heat  capacity  of  the  volatiles  generated  during 
pyrolysis.  Sensitivity  tests  with  the  model  indicated  that  the  shape  of  the 
heat  release  rate  versus  time  curve  was  strongly  dependent  on  the  thermal 
conductivity.  Thermal  conductivity  data  on  this  project  were  limited  to 
temperatures  below  250°C  which  is  the  temperature  at  which  the  lower  molecular 
weight  species  In  wood  start  to  decompose.  It  will  be  necessary  to  extend 
these  data  up  to  700*C  on  char.  The  measurements  must  be  done  in  an  inert 
atmosphere  to  avoid  glowing  combustion  of  the  char.  The  general  formula  for 
the  thermal  conductivity  of  wood,  equation  (4-3),  as  a  function  of  Its  density 
and  moisture  content  is  probably  adequate  at  room  temperature.  However,  the 
temperature  dependence  given  by  equation  (4-6)  does  not  appear  to  be  adequate 
particularly  beyond  its  reported  range  of  validity  which  is  only  up  to  100°C. 
At  the  present  time  it  will  be  necessary  to  measure  the  temperature  dependence 
for  each  species  of  wood  at  its  applicable  density.  It  does  appear,  however, 
that  it  may  be  possible  to  construct  a  temperature  dependent  heat  conduction 
model  in  which  there  is  solid  phase  heat  conduction  through  the  fiber  walls 
and  gas  phase  heat  conduction  and  radiation  through  the  pore  spaces. 

The  general  fornula  for  the  heat  capacity  of  wood,  equation  (4-1),  does 
appear  to  apply  to  all  woods  over  the  temperature  range  for  which  they  are 
stable.  However,  data  on  the  heat  capacity  of  char  at  high  temperatures  is 
needed.  The  heat  capacity  of  the  volatiles  can  be  deduced  from  the  measure¬ 
ments  of  the  heat  capacity  of  the  char. 


Figures  15  and  17  show  che  calculated  heat  release  rates  at  an  external 
flux  of  50  IcW/m^  for  assumed  heats  of  pyrolysis  of  0  and  400  kJ/kg,  respec¬ 
tively.  The  effect  on  the  shape  and  amplitude  of  the  curve  is  significant. 
Therefore,  the  heat  of  pyrolysis  must  be  measured. 

Comparison  of  figures  15  and  20  with  figures  24  and  23  show  that  the 
calculated  ignition  times  are  nearly  twice  as  long  as  the  measured  ones.  It 
should  be  noted  that  the  model  assumes  a  constant  absorptivity  of  0.8  for  the 
Incoming  radiation.  As  the  surface  layer  begins  to  char  the  absorptivity 
should  approach  unity  for  the  real  specimen  in  the  cone  calorimeter  and  thus  a 
shorter  measured  time  to  ignition  would  be  expected.  The  calculated  char 
depth  (see  figure  13)  shows  that  the  calculated  time  at  which  the  surface 
layer  loses  5Z  of  its  mass  is  close  to  the  time  that  the  actual  specimen 
Ignites.  The  calculated  ignition  temperature,  indicated  by  the  cusp  in  figure 
12,  was  420°C  which  is  somewhat  higher  than  the  380*C  reported  by  Quintlere 
and  Harkleroad  [24]  for  Douglas  fir  particle  board.  Choosing  a  critical 
calculated  heat  release  rate  of  less  than  30  kW/a^  would  have  lowered  the 
calculated  ignition  temperature  and  shortened  the  calculated  time  to  ignition. 
As  indicated  earlier,  the  ignition  criterion  in  the  model  must  be  Improved. 

It  should  also  be  noted  that  the  effect  of  diathermacy  is  neglected  in  the 
model.  It  is  difficult  to  judge  its  impact  on  the  time  to  ignition. 

One  of  the  motivating  factors  for  the  introduction  of  a  heat  release  rate 
model  over  that  of  a  mass  loss  rate  model  was  the  variation  in  the  effective 
heat  of  combustion  during  the  burning  period.  However,  the  calculated  effec¬ 
tive  heat  of  combustion  for  dry  Douglas  fir  particle  board  can  be  represented 
by  12  MJ/kg  to  within  *  10Z  over  the  first  three  quarters  of  the  burning 


period  for  all  fluxes  between  25  and  100  kW/m.  The  constancy  of  the  effec¬ 
tive  heat  of  combustion  is  also  demonstrated  by  the  similarity  of  the  calcu¬ 
lated  heat  and  mass  loss  rate  curves  in  figures  14  and  15.  This  is  so  in 
spite  of  a  large  variation  in  the  heat  of  combustion  of  the  vapors  over  the 
full  range  of  decomposition  for  a  single  layer  within  the  specimen.  This 
finding  is  in  agreement  with  the  measurements  in  the  cone  calorimeter.  Thus  a 
constant  effective  heat  of  combustion  might  be  used  in  the  flame  spread 
after  all,  at  least  for  Douglas  fir  particle  board.  However,  it  is  necessary 
to  measure  the  effective  heat  of  combustion  for  the  material  rather  than  use 
the  gross  heat  of  combustion  measured  in  the  oxygen  bomb  calorimeter.  The 
calculations  demonstrate  that  the  low  values  of  the  measured  effective  heat  of 
combustion  of  wood  are  not  due  to  Incomplete  combustion  but  rather  to  the 
difference  in  the  carbon: hydrogen: oxygen  ratio  in  the  vapors  from  that  of  the 
virgin  wood. 

The  formulas  developed  in  Appendices  A,  B,  and  C  are  Included  here 
primarily  because  they  will  be  used  to  determine  the  kinetic  constants  and  the 
heat  of  combustion  of  the  vapors  on  this  project  in  the  future.  The  three 
analyzer  method  used  to  obtain  the  heat  release  rate  and  the  effective  heat  of 
combustion  of  the  specimen  in  the  cone  calorimeter,  as  shown  in  figure  31  and 
32,  serves  as  a  partial  demonstration  of  their  use.  It  is  unfortunate  that 
the  water  vapor  analysis  system  was  not  functioning  properly  during  the 
initial  phase  of  the  project.  The  three  analyzer  method  also  provides  a  means 
of  determining  the  effective  overall  chemical  composition  of  the  volatiles  in 
the  cone  calorimeter  using  the  formulas  of  Appendix  C. 


8.  SUMMARY 


This  is  a  status  report  on  the  prediction  of  the  heat  release  rate  of 
wood.  While  nuch  remains  to  be  done  in  order  to  obtain  a  completely  satisfac¬ 
tory  model,  a  first  cut  on  the  overall  problem  has  been  completed.  A  computer 
program  has  been  written,  laboratory  test  methods  have  been  devised  to  obtain 
the  critical  input  data,  calculations  have  been  made  using  the  model  and 
comparisons  have  been  made  with  experimentally  determined  heat  release  rates 
from-  the  cone  calorimeter.  For  the  limited  comparisons  made,  the  predicted 
and  measured  effective  heats  of  combustion  over  most  of  the  burning  period  are 
within  25Z  and  the  peak  heat  release  rates  are  within  20Z.  The  shape  of  the 
calculated  and  measured  heat  release  rate  curves  are  similar  but  the  time 
scales  are  substantially  different.  This  may  be  due  to  inadequate  data  on  the 
thermal  properties  of  char  at  high  temperatures.  The  second  peak  of  the 
measured  heat  release  rate  curve  is  considerably  smaller  than  the  measured 
one.  This  may  be  due  to  heat  conduction  into  the  marinite  backing  board. 

The  computer  model  is  unique  in  that  it  (1)  accounts  for  the  change  in 
the  heat  of  combustion  of  the  volatiles  generated  during  the  pyrolysis  period, 
(2)  accounts  for  char  shrinkage,  and  (3)  provides  for  different  reaction  rates 
for  the  different  components  of  wood.  However  this  last  capability  was  not 
exercised  during  this  first  test  of  the  model  where  a  single  first  order 
reaction  was  assumed.  This  is  the  standard  assumption  used  in  the  kinetic 
mass  loss  rate  models. 


New  experimental  techniques  have  been  devised  to  measure  the  heat  of 
combustion  of  the  volatile  pyrolysis  products  and  the  thermal  conductivity  of 


wood  and  wood  char  aC  higher  temperatures  than  can  be  attained  for  pyrolyzing 
materials  using  the  standard  methods.  During  the  next  stage  of  the  project 
the  thermal  conductivity  measurements  on  the  char  mist  be  carried  to  higher 
temperatures  still.  Also  high  temperature  data  on  the  heat  capacity  of  char 
must  be  obtained.  The  experimental  data  obtained  on  the  above  properties  were 
determined  from  a  limited  number  of  tests  and  have  not  been  subjected  to  a 
rigorous  error  analysis.  Therefore  they  must  be  regarded  as  preliminary  data 
at  this  point. 

9.  FUTURE  PLANS 

As  this  work  continues  improvements  will  be  made  in  the  model  and  in  the 
experimental  methods  employed  to  obtain  the  input  data.  Furthermore  a  much 
broader  data  base  will  be  obtained  on  the  thermophysical  and  thermochemical 
properties  of  wood  and  wood  char,  particularly  at  high  temperature.  In  order 
to  reliably  describe  the  kinetic  constants  for  the  thermal  decomposition  of 
wood  and  the  effective  heat  of  combustion  of  its  volatile  pyrolysis  products 
as  a  function  of  time,  it  will  be  necessary  to  make  these  determinations  on 
the  individual  wood  components  (i.e. ,  cellulose,  hemi cellulose,  and  lignin). 
The  prediction  techniques  will  also  need  to  be  extended  to  include  other  wood 
products  such  as  fiber  board,  hardboard,  and  plywood.  The  effect  of  fire 
retardants  will  also  be  included. 

Improvements  to  be  made  to  the  computer  model  will  include  (1)  the  reduc¬ 
tion  in  grid  size  and  in  the  time  of  calculation  by  transferring  the  program 
to  a  larger  computer,  (2)  a  state  of  the  art  model  for  the  flame  heat  transfer 
to  the  front  surface,  (3)  an  improved  ignition  criterion,  (4)  modification  of 


the  boundary  condition  at  the  rear  surface  to  simulate  the  marinite  block 
typically  used  to  back  the  specimen  in  the  heat  release  rate  calorimeter  or 
the  simulation  of  the  wall  cavity  during  a  room  fire,  and  (5)  a  model  for 
calculating  the  impact  of  the  fissures. 

In  the  verification  tests  with  the  calorimeter  some  restraint  system  will 
have  to  be  used  to  hold  the  specimen  in  place  as  long  as  possible  in  order  to 
extend  the  period  over  which  the  calculated  and  measured  heat  release  rates 
can  be  compared.  At  present  the  warping  of  the  specimen  near  the  end  of  the 
test  allows  flames  to  move  around  to  the  rear  surface;  or  the  specimen  may 
fall  out  of  the  holder  altogether.  For  selected  tests  the  specimen  will  be 
instrumented  with  thermocouples  to  check  the  calculated  temperature  profiles. 
Some  effort  will  be  devoted  to  measuring  the  heat  flux  from  the  flame  to  the 
surface  in  the  calorimeter. 

Improvements  in  the  pyrolyzer  will  be  aimed  at  better  heat  transfer  to 
the  specimen  to  reduce  the  amount  of  decomposition  that  it  undergoes  by  the 
time  it  reaches  a  constant  pyrolysis  temperature.  Other  problems  with  the 
present  pyrolyzer  will  only  become  apparent  when  it  is  finally  used  in  the  gas 
analysis  mode;  those  tests  should  begin  soon.  It  will  be  beneficial  to  have 
determined  the  critical  thermo chemical  properties  by  two  entirely  different 
methods. 
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Table  1 

Thermal  Conductivity  of  Flberboard  and  Particle  Board 


Thermal  Conductivity  x  10^(lcW/a  .K) 
Flberboard  Particle  Board 


Density 
(kg/m3  ) 


Table  2 


Thermal  Conductivity  of  Oven  Dry 
Southern  Pine  as  a  Function  of  Temperature* 


External 

Heat 

Flux 

(kff/m2) 

Front 

Surface 

Temp. 

<°C) 

Rear 

Surface 

Tfemp. 

CC) 

Average 

Temp. 

Cc) 

Thermal 
Conductivity 
x  104 
(kW/m  .K) 

Rear  Surface 
Convective  Heat 
Transfer  Coefficient 
(kW/mZ.K) 

3.6 

122 

107 

114 

1.28 

.010 

4.6 

139 

121 

130 

1.35 

.011 

5.6 

158 

137 

148 

1.41 

.011 

6.8 

181 

155 

168 

1.36 

.011 

9.1 

213 

180 

197 

1.42 

.012 

12.4 

255 

211 

233 

1.41 

.013 

15.7 

287 

234 

261 

1.46 

0.014 

*Speclmen  Is  1.27  on  thick  and  has  a  density  of  498  kg/m-*. 
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Table  3 


Thermal  Conductivity  of  a  Charred  Specimen  of 
Southern  Pine  as  a  Function  of  Temperature* 


External 

Heat 

Flux 

(kW/m2) 


Front 


Average 


Surface 

Temp. 

(°C) 

Surface 

Temp. 

(°C) 

Temp. 

<8C) 

124 

101 

113 

148 

120 

134 

165 

131 

148 

187 

148 

168 

222 

172 

197 

265 

204 

235 

303 

233 

268 

Thermal 
Conductivity 
x  10** 
(kW/m  .K) 


Rear  Surface 
Convective  Heat 
Transfer  Coefficient 
(kW/m^K) 


is  1.22  mm  thick  and  has  a  density 
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Table  4 


Thermal  Conductivity  of  Oven  Dry  Douglas  Fir  Particle 
Board  as  a  Function  of  Temperature* 


External 

Beat 

Flux 

(ktf/m2) 

Front 

Surface 

Temp. 

ro 

Bear 

Surface 

Temp. 

<#C) 

Average 

Temp. 

(°C) 

Thermal 
Conductivity 
x  10*» 
(kW/m/  .K> 

Bear  Surface 
Convective  Heat 
Transfer  Coefficient 
(kW/m2.K) 

1.8 

83 

74 

79 

1.50 

0.007 

2.7 

109 

96 

103 

1.64 

0.007 

3.6 

133 

116 

125 

1.64 

0.007 

5.4 

170 

148 

159 

1.77 

.006 

7.2 

205 

177 

191 

1.98 

.007 

7.6 

209 

178 

194 

1.85 

.008 

8.9 

233 

199 

216 

1.95 

.007 

11.0 

267 

230 

249 

2.24 

.006 

288 

249 

269 

2.28 

.005 

£  7 

Specimen  .is  1.96  mm  thick  and  has  a  density  of  709  kg/m 
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Table  5 


Thermal  Conductivity  of  a  Charred  Specimen  of  Douglas  Fir 
Particle  Board  as  a  Function  of  Temperature* 


External 

Front 

Sear 

Average 

Thermal 

Bear  Surface 

Heat 

Surf  ace 

Surface 

Temp. 

Conductivity 

Convective  Heat 

Flux 

Temp. 

Temp. 

x  10** 

Transfer  Coefficient 

(kW/m2) 

CO 

<*C) 

(8C) 

(kW/m/  *K) 

(kW/m2.K) 

Table  6 


Input  Data  for 
Property 

Wood 

Activation  Energy 
Pre-Exponential  Factor 
Heat  of  Pyrolysis 
Original  Density 
Heat  Capacity 

Thermal  Conductivity 

Absorptivity 
Contraction  Factors 

Final  Char  Yield 

Volatile  Pyrolysis 
Products 

Heat  Capacity 

Heat  of  Combustion 

Water 

Activation  Energy 
Pre-Exponential  Factor 
Heat  of  Vaporization 
Heat  Capacity 

Flame 

Temperature 

Radiation 

Total  Heat  Flux  to 
Cold  Surface 


Running  the  Model  for  Douglas  Fir  Particle  Board 


Value 


121  kJ/mol 
5.94  x  107  s-1 
0 

709  kg/m3 

1.11  (1  +  0.0067  T)  x 
(0.54  +  0.46Z)  kJ/kg  /.K 

1.24  x  10-4  (0.35  +  0.65Z)  x 
(1  ♦  (T-T0)  6.8  x  10-4)  kW/m/  *K 

0.8  Z  >  0.75;  1.0  Z  <  0.75 

■  Ay  ■  ig  "  1  f°r  o.  65  <  Z  <  1 
-  1  -  (0.65  -  Z)2  for  Z  <  0.65 

0.225 


1.05  +  1.80  x  lO-^T-T^kJ/kg/  -K 
20000  [1.24  -  Z]  kJ/kg 

44  kJ/mol 
4.5  x  103  s”1 
2400  kJ/kg 

4.19  kJ/kg/  .K,  T  <  100 °C 

4.19  +  3.1  x  10“3 
x  (T-100)  kJ/kg/  .K,  T  >  100°C 

1200#C 
10  kW/m2 

30  kW/m2 


Source 


1 

1 

2 

1 

3 

1 

1 

1 

1 


4 

1 

4 

4 

4 

5 


1 

7 

1 


V 


1.  Experimentally  determined  on  project. 

2.  Assumed. 

3.  Temperature  Dependence  -  Koch  [17];  dependence  on  charring  assumed  based 
on  published  value  £or  charcoal. 

4.  Atreya  [13]. 

5.  Handbook  of  Chemistry  and  Physics. 

6*  Calculated. 

7.  Tewarson  [8] 

*  Used  in  basic  computer  program 
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APPENDIX  A 


CALCULATION  OF  THE  DENSITY  OF  THE  GAS  IN  THE  EXHAUST  DUCT 

In  order  to  determine  the  volume  and  mass  flows  in  the  exhaust  duct  using 
an  orifice  plate  or  a  pitot  static  tube  it  is  necessary  to  know  the  density  of 
the  gas.  The  total  mass  flow  through  the  duct  is  given  by 


& 


+ 


(A-l) 


assuming  complete  combustion.  Here  nitrogen  is  taken  to  include  all  of  the 
gases  in  the  atmosphere  except  oxygen,  carbon  dioxide,  and  water.  For 
simplicity  its  molecular  weight  will  still  be  taken  to  be  28. 


Hence  the  density  is  given  by. 

m  _S  „S  „S  .  „S 

p  *S  pN2  xn2  +  0o2  xo2  +  pco2xco2  P^o  xh2o 


(A-2) 


Since 


fS  _  i  _  VS  ^  yS  _  VS 


S  *  1  -  x  o,  -  "'Ho,  - 


equation  (A-2)  can  be  written 


(A-3) 


Hi 


+  (Pn  "  PNJ  Xn .  +  (P 


CO 2  PN2}  XC02  +  ^PH20  ”  XH20 


(  A-4) 


The  density  of  each  species  is  equal  to  its  molecular  weight  divided  by  the 
molar  volume  at  the  reference  temperature.  At  0°C  equation  (A-4)  becomes 


p  -  (28  +  4  Xj  ♦  16  ?  -  10  xj  )/22.4  (A-5) 

A  4  4 

Use  will  now  be  made  of  Che  relationships  between  the  volume  fractions 
measured  in  the  analyzers  and  those  existing  in  the  exhaust  duct.  These  are 
derived  in  Appendix  C.  Since  water  vapor  is  trapped  out  of  the  CO2  analyzer. 


(C-20) 


If  water  vapor  is  trapped  but  CO2  is  not  trapped  out.  of  the  oxygen  analyzer, 


(C-23) 


In  the  cone  calorimeter  both  water  vapor  and  CO2  are  trapped  out  of  the  oxygen 
analyzer  so  that 


xo2  ■  t1  -  42o)  t1  -  4>2)  *S2  <c-25> 

Combining  equations  (A-5),  (C-20)  and  (C-25),  one  finds  that  the  density  of 
the  exhaust  gases  expressed  in  terms  of  the  volume  factions  measured  in  the 
analyzers  is  given  by 

»  -  [4.5  4  [1  -  X^0)  [6.5  -  (1  -  X^)  (4  -  xj  )]1/S.6  CA-6) 

in  the  case  of  the  cone  calorimeter. 


If  the  CO2  is  not  trapped  then  equations  (A-5),  (C-20)  and  (C-23)  are 


combined  Instead  to  yield 


APPENDIX  B 


CALCULATION  OF  THE  FLOW  OF  NITROGEN  THROUGH  THE  EXHAUST  DUCT 


The  formulas  for  Che  calculation  of  carbon,  hydrogen  and  oxygen  release 
rates,  to  be  derived  in  Appendix  C,  require  input  data  on  the  nitrogen  flow  in 
the  exhaust  duct.  The  formulas  for  the  nitrogen  flow  will  be  derived  in  this 
section. 


The  mass  flow  of  gas  through  the  exhaust  duct  of  the  cone  calorimeter  is 
measured  with  an  orifice  plate  using  the  equation 


If  the  orifice  plate  conforms  to  ASTM  specifications,  the  pipe  is  straight  for 

-1  1/2  -1/2 

sufficient  length  on  each  side  and  the  gas  is  air,  C  «  0.031  kg.s  *K  *Pa 
Here  Ap  is  the  pressure  drop  across  the  orifice  plate  in  Pascals  and  Tg  is  the 
absolute  temperature  of  the  exhaust  in  K.  The  orifice  plate  was  calibrated 
(not  on  this  project)  with  a  known  flow  rate  of  CO2  and  C  was  found  to  be 
0.028.  Since  C  is  proportional  to  Che  square  root  of  the  density  of  the  gas 
we  can  define  a  new  calibration  factor, 


C*  •  /p/palr  C 


(B-2) 


and  use  it  in  place  of  C  in  eq.  (B-l).  All  densities  and  volume  flows  in 
these  appendices  are  referred  to  standard  conditions  (760  mm  Hg  and  0°C). 
Nitrogen  is  taken  to  include  all  of  the  gases  in  the  atmosphere  other  than 
oxygen,  carbon  dioxide  and  water.  The  mass  flow  of  nitrogen  is  given  by 


•  -  -  »  •*  -*  -  --  WTJ  «'  .  v-  W_  *_.  *w  v  -.  -.  -»  -.  .  .. 


pN2  \  4/p 


\  _ 

—  X„2  C*  /4P/Te 


Using  equations  (C-20)  and  (C-23)  of  Appendix  C  for  the  case  where  CC>2  is  not 
trapped,  we  find  that  this  can  be  written  in  terms  of  the  analyzer  concentra¬ 
tions  as 


V  •  7===  0  -  4,o)  (i  -  4  -  4 ) c  #  «-4> 

2/pp.  2  2  2  l  e 

K  Mair 

Using  equations  (C-20)  and  (C-25)  of  Appendix  C  for  the  cone  calorimeter  where 
CO2  is  trapped  out  of  the  oxygen  analyzer,  we  find 

v  ■  7=  -  42h*  -  4 c  (B-5) 

2/pp.  2  2  2  J  e 

air 

The  density  of  the  exhaust  gas  is  given  by  either  equation  (A-6)  or  (A-7)  from 
Appendix  A  while  the  densities  p^  and  Pa^r  are  1.25  and  1.29  kg/m 
respectively. 
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APPENDIX  C 


CALCULATION  OF  CARBON,  HYDROGEN,  AND  OXYGEN  RELEASE  RATES  FROM  BURNING  WOOD 

We  wish  to  calculate  the  carbon,  hydrogen  and  oxygen  content  of  the 
volatile  pyrolysis  products  of  wood  from  the  volume  fractions  of  oxygen, 
carbon  dioxide  and  water  vapor  in  the  combustion  products  along  with  the  mass 
flow  of  nitrogen  through  the  system.  It  is  assumed  that 

1.  Wood  is  composed  only  of  carbon,  hydrogen,  and  oxygen. 

2.  All  of  the  volatiles  undergo  complete  combustion  to  water  vapor  and 
carbon  dioxide. 


There  are  two  situations  of  concern.  The  first  is  the  pyrolyzer  where 
the  pyrolysis  takes  place  in  the  nitrogen  stream  and  then  is  swept  into  a 
methane  flame  where  complete  combustion  takes  place.  The  second  situation  is 
the  cone  calorimeter  where  there  is  no  supply  of  methane  and  the  volatiles 
simply  burn  in  their  own  flame.  The  calculational  procedures  are  identical. 
The  flow  rate  of  methane  is  simply  set  equal  to  zero  for  the  cone  calorimeter 
but  is  measured  for  the  pyrolyzer  or  simply  calculated  from  the  volume 
fractions  of  O2,  CO2  and  ^0  in  the  exhaust  when  the  wood  is  absent. 


The  water  vapor  volume  fraction,  Xg  Q,  is  the  same  in  the  analyzer  as  it 


is  in  the  exhaust  duct.  However,  the  O2  and  CO2  volume  fractions. 


X^  and  ,  in  the  analyzers  differ  from  their  values  in  the  exhaust 

*  s  L  s 

duct,  Xn  and  Xrn  ,  due  to  trapping  of  water  vapor  from  the  analyzers. 

2  .  C02 

X^ 


Furthermore  Xq  will  depend  on  whether  COj  is  also  trapped  out  of  the  oxygen 


analyzer  as  it  currently  is  in  the  cone  calorimeter.  Since  both  cases  arise, 
the  following  calculations  will  consider  Case  I  where  CO2  is  not  trapped  and 
Case  II  where  it  is  trapped.  It  will  be  assumed  that  the  molar  flow  ratios 
are  the  same  in  the  analyzers  as  they  are  in  the  exhaust  duct  for  those 
species  which  are  present  in  both. 


The  mass  flow  rates  will  first  be  calculated  in  terms  of  the  volume 
fractions  in  the  stack  or  exhaust  duct.  Then  they  will  be  transformed  for 
each  of  the  two  cases  to  be  expressed  in  terms  of  the  volume  fractions  in  the 
analyzers.  The  stack  volume  fractions  are  given  in  terms  of  the  molar  flows 
by 


C0„ 


4,0 


S' 

•  *  • 

*-  +  n_-  - 

°2  C02 

hV 

°co2 

CM 

JP 

1 

+  n  +  & 

o2  co2 

+  ^h2o 

“h2o 

4 

+  f»+  hnn 
(>2  CO2 

+  ^h2o 

4 

and 


2  °n2  +  no2  +  nco2  +  &h2c 


(C-l) 


(C-2) 


(C-3) 


(C-4) 


Furthermore 


Here  nitrogen  Is  taken  to  Include  all  of  the  gases  in  the  atmosphere  other 
than  oxygen,  carbon  dioxide,  and  water.  For  simplicity  its  molecular  weight 
will  still  be  taken  to  be  28. 


The  molar  flows  of  the  combustion  products  in  the  duct  are  related  to  the 
carbon,  hydrogen,  and  oxygen  content  of  the  volatile  pyrolysis  products,  the 
species  present  in  the  incoming  air  and  the  methane  by  the  expressions. 


nC02  "  nC  +  nC02  +  nCH4 


(06) 


.  1.  .  .  O  A  . 

^0  “  2  ^  ^0  2  °CH4 


(07) 


n02  “  2  &0  +  &02  fhC02  ~  &C02^  “  2 


(08) 


Combining  equations  (C-l)  to  (C-8)  the  molar  flow  ratios  of  atomic  carbon, 
hydrogen,  and  oxygen  in  the  volatile  pyrolysis  products  to  the  molecular 


nitrogen  flowing  through  the  system  are  given  by 
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The  ambient  volume  fractions  x£  ,  x£Q  ,  and  X^  ^ 
analyzers  prior  to  the  test  with  the  traps  bypassed. 


are  determined  from  the 


Noting  that 


3  nC 

\ 

"  7  . 

% 

.  1  *s 

s 

"  IS  . 

*0 

4  n0 

X 

“  7  . 

X 

and 


4^4 

• 

“  7  . 

\ 

nN2 

(C-12) 


(C-13) 


(C-14) 


(C-15) 


we  find  that 


(C-16) 


"H 


and 


®N2"  4  aCH4 


(c-i/; 


“o 


4 

7 


2 4>„+ 2 4* 4,o  2 4  ♦  2 4,+ 


1  ~  -  X? 

*o2  aco2  ^0 


1  —  X°  —  X°  —  X? 

°2  C02  V 


V, 


(018) 


Equations  (016),  (017)  and  (C-18)  are  in  terms  of  the  volume  fractions  in 
the  exhaust  ducts.  It  is  necessary  to  establish  the  relationships  between  the 
volume  fractions  in  the  analyzer  and  in  the  exhaust  duct. 


The  concentration  in  the  C02  analyzer  is  given  by 

nco 

-  - - - -  (019) 

2  nN2  n02  nC02 

where  the  molar  ratios  are  taken  to  be  the  same  in  the  analyzer  as  they  are  in 
the  exhaust  duct  for  N2,  02,  and  C02«  Combining  equations  (02)  and  (019) 
yields 


(*  -  42o)  42 


(020) 


If  C02  is  not  trapped  out  of  the  oxygen  analyzer  (Case  I)  then 


If  CO2  Is  trapped  (Case  II)  then 


A'  1 


(C-22) 


Combining  equations  (1)  and  (21)  yields 


A1  -  (» -  42o)  42 


(C-23) 


for  Case  I.  Combining  equations  (C-l)  and  (C-22)  yields 


xo2  ■  (*  -42  -  42o)  *5’ 


(C-24) 


for  Case  II.  Combining  equations  (C-20)  and  (C-24)  expresses  the  volume 
fraction  of  02  in  the  exhaust  duct  in  terms  of  the  volume  fractions  in  the 
analyzers. 


*02  -  t1  -  42o)(‘  -  4o2)  ’v 


(C-25) 


It  can  be  seen  from  equations  (C-20)  and  (C-23)  that  for  Case  I 


1  "  X02"  XC02~  ■  C1  “  *H20^1  “  *02" 


(C-26) 


Substituting  (C-20),  (C-23)  and  (C-26)  into  (C-16),  (C-I7)  and  (C-18)  yields 
for  Case  I  with  no  trapping 


^o  -xco 


"  1  1-X°  -x° 

coj  0,  co„ 

2  2  2 


yO  1%  A  “c 


.  J  xX°  2{^1]  WA .  (c_29) 

■°  7  U-42o)f-42-^o2)  (K  -^)  -O^O^O )  ^ 


It  can  be  seen  from  equation  (C-20)  and  (C-25)  that 


1  '  ^  *  42  '  ^2o  *  (‘-^oX1-^1-^! 


(C-30) 


Substituting  equations  (C-20),  (C-25)  and  (C-23)  into  equations  (C-16),  (C-17) 
and  (C-18)  yields  for  Case  II  (with  trapping) 


.  3 


v'7  <c'31) 


.  3  . 


(  1-XH20^  ( ( 1_XC02^ 


2  X-  +  2  X- 


(i-x£  )(i-x£0  )  (i-4.0)(i-xj JCi-xio  ) 


Now  we  drop  the  superscripts  S,  A  and  A'  and  assume  that  the  volume  fractions 


are  those  measured  in  the  analyzers,  (X^  Q,  is  the  volume  fraction  in  the 


water  vapor  analyzer).  We  are  assuming  Case  II  where  the  CO2  is  trapped. 
Then 


X, 


.  3 


CO* 


X 


'CO* 


u  -  «* 

C  7 


2 _ \ .  3  . 

,0  „o  „o  IV  4  mCH, 


(C-34) 


S'-hJV-hoJ  1-xo2-xco2-xh2o /  2  "  w“4 


(1-XH20,tl'X<>2)(1'XC‘>2) 


1-xo2-xoo2-xS2o, 


V  i  V  (c-35) 


4 1 2  (V  V  , 

“o  "  7  ' - + 


V 


2(XOj<'XC02)'fXH2o' 


(1-xo2Ki-xoo2)  (‘-xh2oK1-xo2)(1-xco2)  I  "2 
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•o  Method  of  section  4.3 
□  Guarded  hot  plate 
&  Calculated  from  Eq.  (4-6) 
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Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Pir  Particle  Board  Exposed  at  an  External  Radiant 
Flux  of  25  kW/m2 


Measured  Effective  Heat  of  Combustion  for  Dry  12.7  mm  Specimen 
of  Douglas  Fir  Particle  Board  Exposed  at  an  External  Radiant 
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Measured  Heat  Release  Rate  Deternined  by  the  Three  Analyzer 
Method  for  Dry  12.7  mn  Specimen  of  Douglas  Fir  Particle 
Board  Exposed  at  an  External  Radiant  Flux  of  50  kW/m2 
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